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In this paper we study the main geometric properties of the Carnot-Caratheodory 
(abbreviated CC) distance dn in the setting of fc-step sub-Riemannian Carnot 
groups from many different points of view; see Section 12.11 and Section 12.21 An 
extensive study of the so-called normal CC-geodesics is given. We state and prove 
some related variational formulae and we find suitable Jacobi-type equations for 
normal CC-geodesics; see Section 12.21 One of our main results is a sub-Riemannian 
version of the Gauss Lemma; see Section 12.41 We show the existence of the metric 
normal for smooth non-characteristic hypersurfaces; see Corollarv I2.31[ In Section 
12.31 we compute the sub-Riemannian exponential map exp ^ for the case of 2-step 
Carnot groups. Other features of normal CC-geodesics are then studied. In Section 
12.61 we show how the system of normal CC-geodesic equations can be integrated 
step by step. Finally, in Section [3] we show a regularity property of the CC-distance 
function 5h from a C'^-smooth hypersurface S; see Theorem 13.51 
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1 Introduction 



In the last thirty years, many progresses have been done in developing Analysis through 
very general geometries and metric spaces. This trend in mathematical research, already 
initiated by Federer's treatise [21], has been pursued by many authors, with different point 
of views: Ambrosio [3] , Ambrosio and Kirchheim [H [5] , Capogna, Danielli and Garofalo 
[To], Cheeger [12], Cheeger and Kleiner [13], David and Semmes [18], De Giorgi [19] . 
Gromov [35l[36], Franchi, Gallot and Wheeden [25], Franchi and Lanconelli [26], Franchi, 
Serapioni and Serra Cassano [271 130], Garofalo and Nhieu [31], Heinonen and Koskela [37] , 
Jerison [30j, Korany and Riemann [l2], Pansu [56] but, of course, this list is not complete. 

Sub-Riemannian or Carnot- Caratheodory geometries have become a research field of 
great interest, also because of their wide connections with different fields in Mathematics 
and Physics, such as PDE's, Control Theory, Mechanics, Theoretical Computer Science. 

For references, comments and some different perspectives on sub-Riemannian geometry 
we refer the reader to Agrachev and Ghautier [1], Agrachev and Sarychev [2], Gromov 
[36] , Montgomery [54j , Pansu [561 E3 ; Strichartz [61] , Vershik and Gershkovich [62] . 

Very recently, the so-called Visual Geometry has also received new impulses from this 
field; see [H], [15] and references therein. 

In this article we begin a study of fine properties of the function "sub-Riemannian 
distance from a hypersurface" in Carnot groups. Before giving a formal outline of the 
context and content of the paper, we would like to briefly and colloquially discuss what 
we think are the main features of this research. 

The elementary analysis of the function "distance from a closed set" E with smooth 
boundary dE in the Euclidean space is based on the intuitive idea of "moving" a 
sufficiently small closed Euclidean ball BmW^\E until it touches dE at some point P. 
Let -Bo be the ball B in its final position and let O be its center. The Euclidean segment 
OP is (part of) the segment perpendicular to dE at P. For points Q in OP, the distance 
from Q to is realized by the segment OQ. The idea of a "sufficiently small moving ball", 
however intuitive, is not easy to deal with in calculations and it is not really necessary. 
In fact, it is easier to work with the final configuration: given a point P on dE and the 
direction u normal to dE at P, we can recover Bq as any of the sufficiently small balls 
centered at a point O on the segment leaving P in the direction i/, having radius OP. 
The normal direction z/ is a first order differential object; while the fact that the ball Bq 
touches E at P only, generally depends on the "extrinsic" curvatures of dE at P, which 
is a second order differential object. 
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Exactly the same line of reasoning works if E lives in a Riemannian manifold. The 
direction z/, normal to dE at P, determines the geodesic 7 having speed u at P, the 
geodesic normal to dE at P. The ball Bq can then be recovered as a ball having center 
on 7, sufficiently close to P. 

The sub-Riemannian case presents new difficulties. Consider, for expository reasons, 
the case of M"^ endowed with a sub-Riemannian structure of codimension f > 1. In 
particular, to each point x G is associated a linear space C T^M."" of admissible 
horizontal directions, dim(iJj.) := h = n — v, and a metric to measure lengths of 
vectors in H^. We stress that the metric Qh can also be defined (non-canonically) as the 
restriction of a Riemannian metric g defined on the whole tangent bundle TM". 

We can use this structure to define lengths of admissible curves (those whose velocity 
belong to H) and, under suitable hypotheses, this leads to a new distance in M", 
called Carnot-Caratheodory (abbreviated CC) distance. This distance du is a length- 
metric, realized by geodesies In this introduction, we assume they are smooth, although 
the regularity of the so-called CC-geodesics is, in general, an open problem; see Section 
Oand Section EJl 

We can mentally repeat the previous procedure, moving a small metric ball B until it 
touches dE at x. We assume that x is non- characteristic, i.e. roughly speaking, the linear 
space Hx is not totally contained in T^dE, the tangent space to dE at x. In particular, 
we can determine a direction i^^ (x) in (the unit horizontal normal to dE at x) which 
is orthogonal with respect to gn^ to all directions in fl T^dE. The novelty is that there 
are infinitely many CC-geodesics 7 passing through x and having tangent vector z/^(x) 
there. In fact, there is a f-parameter family of them, where v = codimif. Hence, we can 
not identify the final ball Bq based on the information contained in alone. 

However, in several cases there exists a unique metric normal 'fj^j- to dE at x. By this 
we mean the (maximal) CC-geodesic curve 7_v' starting from x, with the property that, 
for each y in •y^ \ E, dH{x,y) realizes the distance between y and E. i.e., the metric 
ball centered at y and having radius du {x, y) touches dE at x. The notion of metric 
normal was introduced by the first two authors in the context of the (lowest dimensional) 
Heisenberg group H^; see [H]. 

In order to specify the metric normal among the CC-geodesics tangent to aX x & dE 
we need i real parameters. In the Heisenberg case, « = 1, the extra parameter could be 
read off from the parametrization of dE. This fact is crucial, because from the equation 
for dE one can explicitly write down the "exponential map" associated with dE, and, for 
instance, deduce regularity properties of the function "distance from . 

In [7J it was observed that, in the Heisenberg case, the extra parameter could be read 
as a "imaginary curvature" for dE. At the same time and independently, the third author 
|51j . continuing his work [50j on variational formulas in sub-Riemannian Carnot groups, 
associated to each non-characteristic point x G dE a "vertical vector" w of parameters, 
which reduces to the imaginary curvature in the Heisenberg case. These parameters, 
which can be interpreted as "curvatures" (or Lagrangian multipliers; see Section 12. 2p . 
have no immediate counterpart in Riemannian geometry and play an important role in 
the geometric analysis of dE in this sub-Riemannian setting. An interesting feature of 
these "curvatures", which appears in seminal form in [6] and in much greater generality 
in |50[ [51], is they can best be computed in terms of the Riemannian gradient (with 
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respect to the fixed Riemannian metric g on the base manifold) of the smooth function / 
defining dE locally (i.e., dE = {/ = 0}, locally). Roughly speaking w can be regarded 
as the "non-horizontal" (more precisely, vertical; see Section [LTI) part of the Riemannian 
gradient of /; see Definition II. 301 

Researchers working with "first order" geometric measure theory of hypersurfaces in 
Carnot group (isoperimetric inequalities, Lipschitz submanifolds, functions of //-bounded 
variation, etc.) are generally concerned with the horizontal gradient of the given defining 
function only. 

It is in fact an accepted rule-of-thumb that horizontal vectors contain all the first order 
geometric information of the group, and that this extends to some second order objects 
(e.g., //-mean curvature and minimal surfaces). When dealing with such "second order 
information" as that encoded by the distance function, the "forgotten component" of the 
Riemannian gradient has to be taken into account. In this paper, we bring together these 
two different perspectives, relating these "curvatures" w to the metric normal 7_^, hence 
to the distance function from a smooth hypersurface. 

This way, among others things, we shall extend to Carnot groups many results of [6] 
and [7], but with a somewhat different approach. 

We now give a less informal view of the paper and of its context. 

Definition 1.1 (Metric Normal). Let {X,d) be a metric space and let E be closed in X. 
The metric normal to E at x & E is the set := {y E X : d{x,y) = d{y,E)}, where 
d{y,E) := mi^(zE d{z,y) . 

In this paper, the metric space {X, d) is a Carnot group G endowed with its "natural" 
CC-distance dn ■ 

A k-step Carnot group G is a Lie group, whose Lie algebra g admits a stratification 
Q = H ® V, where V := ©^^=2^fc — H, H^+i = [Hj^H^] and i/fc+i = {0}. In H we 
find the so-called horizontal vectors, while in V we find the vertical vectors. 

Notation 1.2 (Projections). Throughout this paper, the mappings Vh^ : TG — > Hi, 

Vv '■ TG — > V , will denote the projection operator onto the subbundles Hi {i = 1, ...,k) 
and V , respectively. 

Remind that, as for any Lie group, an internal group of translations, called "left- 
translations", is given on G. In addiction, any Carnot group G has intrinsic dilations 
making it a homogeneous group; see Section [LTI In general, we shall say that a property 
in G is intrinsic whenever it is invariant with respect to these transformations. 

The horizontal space H is endowed with a Riemannian metric gn which is used in the 
obvious way to define lengths of horizontal curves 7 in G (i.e. 7 : M — > G is absolutely 
continuous and 7(t) G H for a.e. t G M). The distance between points in G is the infimum 
of the lengths of curves joining them, and it can be proved that such distance is realized 
by the length of (not necessarily unique) CC-geodesic curves; see Section 12.11 and Section 
12.21 for a precise definition of CC-geodesic. 

We stress that is a difficult open problem showing, or disproving, that CC-geodesics 
are smooth for all Carnot groups; see ^54^ i55j. 
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In the Euclidean case, by assuming reasonable hypotheses on dE, the metric normal 
at a point x G dE can be identified with the normal direction associated with the surface 
at X G dE, because it is a segment; see [2T], [22]. As a consequence, whenever we recall 
the normal direction, we emphasize the linear aspect of the metric normal, forgetting 
that it is, first of all, a path. In the Riemannian setting it is known that, if we stay 
near to a smooth hypersurface 5* there exists one, and only one, unit vector belonging to 
the tangent space to the base Riemannian manifold that select the metric normal to the 
hypersurface S at x. As already said, in the sub-Riemannian setting this problem it is 
not so straightforward. Indeed, just to fix the ideas in the simplest Heisenberg group 
(see ^52j for further details), even when we fix a non-trivial horizontal vector z/^(x) G Hx, 
wher^ H = spauj^jX, Y}, infinite different CC-geodesics starting from a; G exist with 
the same initial velocity (x). In [6j this problem was solved by remarking that a path that 
locally parameterizes the metric normal of a smooth surface S at x G S" (whenever X is 
non-characteristic) can be selected just by considering two intrinsic object associated with 
the surface S at x, i.e. the intrinsic unit normal vector z/^(x) and the so-called "imaginary 
curvature" of S at x; see |6]. Roughly saying, if S* = {y G : f{y) = 0}, with / G C^{M^) 
and X = [xi,yi,ti] G S* is non-characteristic, i.e. Xf{x) = dx^f\x) + 2yidt^f{x) 7^ or 
Yf{x) = dyj{x) - 2xidtJ{x) ^ 0, where x = [xi,yi,ti], then 



(which are, respectively, the intrinsic unit normal along S* at x and the imaginary curvature 
of S at x) select an horizontal path that is a subset of the metric normal jj^f to S at x. 

The problem we trait in this paper concerns the generalization of this program to any 
Carnot group. Later on we shall describe our approach. 

Assume that S is the boundary of an open bounded subset f2 of a Carnot group G. 
Now let us assume that the set Q satisfies the so-called property of the internal ball. 
Namely, any point x G 5" can be touched by a closed metric ball B{y,r) C fl and 
dB{y,r) H S = {x}. In this case dB{y,r) and S have the same tangent space at x. The 
objects we are searching for, whenever they exist, are somehow "hidden" in the invariants 
associated with this tangent space. On the other hand the CC-ball B{y, r) is the union of 
all the horizontal paths of length less than r, starting from y. This fiux of CC-geodesics 
is usually determined by the solutions of a suitable Hamiltonian system; see Section 12. 1[ 

At this point we have to select the path 7 : [0, r] — > G, 7(0) = y, 7(r) = x, connecting 
y to X and satisfying the property dn (x, 7(t)) = mf;,^QB(y,r) dn (7(^)5 z). 

We will need two key facts: (i) the CC-distance from a point satisfies the eikonal 
equation (see (5^), i.e. Igradifdnl = 1, at each regular point of dn] (ii) even in this 
sub-Riemannian structure, a Gauss-type lemma holds (see Section [274|) . namely if 7 is a 
normal CC-geodesic leaving the center of the ball, then ■j(t) = gradndHijit)). We stress 
that the last identity is a straightforward consequence of the eikonal equation. 

The remaining part of the proof can be described as follows. 

^According to the notation used in [51 [7], in this introduction wc adopt the foUowing convention: every 
point X G is given, using exponential coordinates, by the triple x = [xi,yi,ti]. Moreover, we use, 
as a vector basis {X, F} for the horizontal space H, the following vectors field: X{x) := d^-^ + 2yidti, 
Y{x) :— dy^ — 2xidti- This convention for will be slightly changed in the sequel. 



X/(x)X(x) + r/(x)F(x) 
v/(X/(x))2 + (F/RF 



v/(X/(x))2 + (F/(x))2' 



[X,Y]f{x) 
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Let us introduce a left-invariant frame X_ := {Xi, . . . , X/j, Xh+i-, • • • , X„} for TG where 
H = spanjgjXi, . . . , X^} and F = spauj^jX/j+i, . . . , X„}. Moreover we fix a Riemannian 
metric ■) := (■, ■) on TG that makes orthonormal the frame X. We shall assume that 
the restriction oi g to H equals the sub-Riemannnian metric Qh, i.e. gin = Qh ■ 

Set Piit) ■= {XidH){n{t)) {I = l,...,n), where dni-) := d„{-,y). 

In particular, one has 

Pi{t) = {grad{XidH){l{t)),^{t)) {I = 
By the Gauss' lemma (i.e. 7(t) = gradjjdH {^{t))), we get 

h 

m = Y,{x,xi){dnhmpAt)- 

On the other hand X^X, = [X^-, X;] + X^X^, thus 

h h 

3=1 i=i 
Here plays a role the eikonal equation, namely X]j=i(^i'^«)^(7(^)) — 1' because 

m = ^X, (X^(XA)'(7W)P,W)+E[^^'^'](^-(^W))^:'W 

Vi=i / i=i 

h 

= Y.[x„xi]{dn{impAt)- 

It turns out that [Xj, Xi]{dH {'jit))) = X]a=/t+i ^j't-^ai^H^'jit))), where the coefficients 
are the structural constants of the Lie algebra g; see Section II. 1[ 
This way, we obtain the following system of O.D.E.'s: 

h n 

^'W = E E CfiP,{t)P^{t) {l = l,...,n). (1) 

j=l a=h+l 

By using the properties of the Carnot structural constants (in particular, see ([6])), this 
system can be solved step by step. A similar algorithm will be discussed in the Appendix 
of Section [2J 

The previous discussion, for which we refer the reader to Section [231 yields the next: 

Lemma 1.3. Let 7 : [0,r] — > G (r > 0) be any normal CC-geodesic^ of unit-speed 
and parameterized by arc-length. Let 7(0) = y, P{0) = Pq be its initial data and set 
dn (x) = dn {x, y) (x G G) . Then we have 

(i) gradndHi'jit)) = Pnit) for every t G [0,r]; 

^See Section l2.ll for a precise definition of normal CC-geodesic. 
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(a) gradydn {'y(t)) = Pv{t) for every t G [0,r]; 
(m) m = Y.U Tl=H+i CI P,it)Pait) for every t e [0, r] 



{l = l,...,n). 



This lemma solves the problem of selecting, for any regular point x belonging to the 
CC-sphere r) := dB{y,r), the unique normal CC-geodesic having velocity vector 

equal to the horizontal normal direction at that point (i.e. Ph (0) = ^^(a;)) and connecting 
this point to the center y of the CC-sphere. In fact, by uniqueness of solutions of O.D.E.'s 
systems, one gets that the desired curve must be the normal CC-geodesic defined by 



Here we have set 



where u denotes the Riemannian unit normal along the CC-sphere §^(?/,r), thus is 
just the horizontal unit normal along §^(|/, r) and zu = ^j^- Moreover, exp^^ denote 
the sub- Riemannian exponential map; see Section 12.31 
An immediate "geometric" consequence can be given: 

Corollary 1.4. Let S = {x & G : f{x) = 0}, where f is a C'^ function. Assume that 
there exists a CC-ball B{y,r) with center y and radius r, such that B{y,r) C {f{x) < 0}, 
or B{y, r) C {/(x) > 0}, and B{y, r) fl 5 = {x}, where x E S is non- characteristic. Then 
there exists the metric normal 7_^^ to S at x and for every t G [0,r], 7(t) G 7^^, where 



Till now we have described some of the main results of this paper, but there are many 
others aspects and related questions. Here we briefly give an account of the rest of the 
paper. 

In Section 11.11 we shall introduce the main ingredients to deal with Carnot groups 
and their sub-Riemannian structures: Lie algebraic preliminaries, Carnot dilations, CC- 
distance, properties of the Carnot structural constants and the notion of curvature of a 
distribution. Moreover we shall discuss some other tools as, for instance, the Levi-Civita 
connection, the so-called iZ-connection (or horizontal connection and, the related notion 
of i?-connection) and covariant differentiation along curves. We also give some basic 
examples. 

In Section 11.21 we will just recall notation and some basic definitions to work with 
hypersurfaces of Carnot groups, such as the notion of characteristic point and that of 
horizontal perimeter. 

In Section [2] we begin our study of Carnot-Caratheodory metrics and their geodesies. 

In Section [2m we will introduce CC-geodesics from the so-called Hamiltonian point of 
view, well-known in literature; see [8], [51], [31]. We will give the equations for normal 
and abnormal curves (and minimizers) and discuss some examples. Here we just remind 
the so-called normal CC-geodesic equations: 



7(t) := exps^iy, -Afix)){t) 



t G [0,r]. 



7(t) := exps^{y,-Af{x)){t) 



t G [0,r]. 



(Normal Equations) 
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where P = {Ph ,Pv) is the n- vector of the momentum functions associated with a fixed 
orthonormal (left-invariant) moving frame X. = {^i,---jXn} for TG. Furthermore, 
C{Pv) G Ain is a hnear combination of n x constant matrices which only depends 
on the Carnot structural constants. More precisely, one has 

n 

C{Pv) := J2 

a=h+l 

where := [C^j]ij=i^,,,^n {a = h + 1, ...,n) and, by definition, C°j := {[Xi,Xj],Xa). 

In Section 12.21 we shall discuss another "natural" point of view on this subject: the 
Lagrangian one. This is because it seems the more natural way to obtain some additional 
information about normal CC-geodesics and their minimizing features. In particular, 
we will derive both the first and the second variation of the natural sub-Riemannian 
Lagrangian. For precise statements see Proposition 12.131 and Theorem 12.141 Then, after 
introducing the notion of CC-geodesic variation, we shall prove the validity of another 
interesting second variation formula; see Corollary 12.171 Starting from these formulae, 
we will deduce the natural Jacobi equations for normal CC-geodesics; see Definition 12.181 
Furthermore, we will discuss another Jacobi-type system of O.D.E.'s, which is obtained 
under suitable assumptions on the variations. These "restricted" Jacobi-type equations 
read as follows: 

Vf y + R{Ph , Y)P„ + C{Pv )VtY = 0, 

where Vt denotes "covariant differentiation" and R denotes the Riemannian curvature 
tensor. This material can be useful in the study of the conjugate and cut loci of a point, 
in this context. 

In Section [ZS] we shall define the sub-Riemannian exponential map exp^ and we will 
show some of their basic features. 

In Section 12.41 we will prove a sub-Riemannian version of the Gauss Lemma; see 
Proposition 12.291 and Corollary 12.311 

In Section 12.51 we will compute the sub-Riemannian exponential map exp ^ for the 
special case of 2-step Carnot groups. Remind that the sub-Riemannian exponential map 
based at the point Xq G G is a mapping 

exp5^(a;o,-)(-) : UH x H2XR — ^ G, 

where UH = {X G H : |X| = 1} denotes the bundle of all horizontal unit vectors. We 
will show that 

exps„{xo,Po){t) ■.= xo+ [ e-^^^^^^^^^PH^ds- ^ V f / {C^Xn ,Xh) ds^ e^^, 

^ a=h+l V -^0 / 

where Pq = {Ph{<^),Ph,) e UH x H2, C% := [C^]i,j=i,...,h E Mh {a = h + l,...,n) and 
CniPH,) := EL^+i^^C*-. Moreover 

XH{t):=XH{0)+ [ e-^«(^«2)^p^(o)rfs. 
Jo 
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In the previous formula we have used the notation for the exponential of a linear operator. 
More precisely, e"*"^ ^^"2 denotes the exponential of a square-matrix, i.e. the hxh-matrix 
defined by: 

e-» '".^ >■ :^ Id. - C, (P,. ), + l^^ifiW! _ lC.(P.)sf ^ 

We shall analyze the existence of T-periodic solutions for the "auxiliary" horizontal path 
Xh (t) = {xi{t), ...,Xhit)) previously defined. Indeed, the T-periodicity of Xh (t) is somehow 
connected with the study of the conjugate and cut loci of a point; see Remark [2.411 

In Section [2l6] we will show how, at least in principle,the system of normal CC-geodesic 
equations can be integrated step by step. 

In Section [3] we will apply some of the tools previously developed toward the study of 
the CC-distance function 6h from a C'^'-smooth {k > 2) hypersurface S, i.e. 

6h{x) := inf dH{x,y). 

This will be done only for 2-step Carnot groups, by using the explicit structure of exp <^ 
in this case. More precisely, we shall define a mapping $ : Sx] — e, e[ > G (e > 0) by 

^y,t) ■.= exps^{y, ^iy)){t), 

and then we will compute its Jacobian; see Lemma 13. 2[ To be more precise, let us fix 
Riemannian normal coordinates [ui, u„_i) around yo G 5*. Then V{y) := ^ A ... A ^ 
is a normal (non-unit) vector along S, in a neighborhood oi yo & S and it turns out that 

|det [Jiy,o)^] \ = \VhV\, 

where x7{yfl)^ denotes the Jacobian matrix operator at (y, 0) G Sx] — e, e[. Therefore, out 
of the characteristic set Cs, the map $ turns out to be invertible. 
Finally, we shall prove the following (see Theorem 13.51) : 

Theorem 1.5. Let G be a 2-step Carnot group. Let S G G be a C^-smooth hypersurface 
with k > 2 and let 6h denote the CC-distance function for S . Set Sq := S\ Cs, where Cs 
denote the characteristic set of S . Then, for every open set Uq compactly contained in Sq, 
there exists a neighborhood U G G of Uq having the unique nearest point property, 
with respect to the CC-distance. Furthermore, the CC-distance function from Uq (1 S is 

\u\Uo 



^h\u\Uo -smooth function. 



The proof of this result is based on explicit computations and on the sub-Riemannian 
Gauss Lemma (see Proposition I2.29( ). 



1.1 Sub-Riemannian Geometry of Carnot groups 

In this section, we will introduce the definitions and the main features concerning the 
sub-Riemannian geometry of Carnot groups. References for this subject are, for instance, 
do], [32], [31], [351 [36], [13116], [19], [M], [56l[57], First, let us consider a C°°-smooth 
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connected n-dimensional manifold and let H C TN be a /ii-dimensional smooth sub- 
bundle of TN . For any p & N, let denote the vector subspace of TpN spanned 
by a local basis of smooth vector fields Xi{p), ...,Xhi{p) for H around p, together with 
all commutators of these vector fields of order < k. The subbundle H is called generic 
if for all p G dimT^ is independent of the point p and horizontal if = TN for 
some G N. The pair (A^, H) is a k-step CC-space if is generic and horizontal and if 
k := inf{r : T^ = TN}. In this case, we have that 

= T^ C H = T^ C T^ C ... C T'' = TN (2) 

is a strictly increasing filtration of smooth subbundles of constant dimensions rii := dimT' 
{i = l,...,k). Setting {Hi)p := T^ \ T'-^, then gr{TpN) := ®\^^{Hk)p is the associated 
graded Lie algebra, at the point p E N, with Lie product induced by [■,■]. Moreover, 
we shall set hi := dimi/j = — ni_i (no = ho = 0). The A;-vector {hi,...,hk) is the 
growth vector of H. Notice that every Hi is a smooth subbundle of the tangent bundle 
TT : TN — > N, i.e. vr^ : Hi — > N, where vr^. = 7r\Hi {i = 1, k). 

Definition 1.6. We will call graded frame X. = {^i, ^n} for N, any frame for 
N such that, for any p E N we have that {Xj.(p) : uj^i < ij < Uj}, is a basis for 
H,\p:=H,{p) (j = l,...,fc). 

Definition 1.7. A sub-Riemannian metric Qh = {■, ■)h on N is a symmetric positive 
bilinear form on H. If {N,H) is a CC-space, then the CC-distance dH{p,q) between 
p, q E N is 

dH{p,q):=mf j a/ {'j,'j)„dt, 
where the infimum is taken over all piecewise- smooth horizontal paths 7 joining p to q. 

In fact. Chow's Theorem (see [35], [51]) implies that dn is actually a metric on A^, since 
any two points can be joined with (at least one) horizontal path; moreover the topology 
induced by the CC-metric turns out to be compatible with the given topology of A^. 

The general setting introduced above is the starting point of sub-Riemannian geometry. 
A nice and very large class of examples of these geometries is represented by Carnot groups 
which for many reasons play, in sub-Riemannian geometry, an analogous role to that of 
Euclidean spaces in Riemannian geometry. Below we will introduce their main features. 
For an introduction to the following topics, we suggest Helgason's book, [38], and the 
survey paper by Milnor, [17], regarding the geometry of Lie groups, and Gromov, [35] . 
Pansu, [551IFF] . and Montgomery, [51], specifically for sub-Riemannian geometry. 

By definition a k-step Carnot group (G, •) is a n-dimensional, connected, simply con- 
nected, nilpotent and stratified Lie group (with respect to the group law •). This means 
that its Lie algebra g = M" satisfies: 

g = Hi®...®Hk, [Hi,H,^i]=H, {t = 2,...,k), Hk+i = {0}, (3) 

We shall denote by the identity on G. Any x G G defines smooth maps Lxj Fix '• G — >■ G, 
called left-translation and right-translation, respectively, by Lx{y) := x»y, Rx{y) '■= yx, 
for any y E G. Remind that the Lie algebra g is naturally isomorphic to TqG by identifying 
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any left-invariant vector field X with its value at 0. The isomorphism is explicitly given 
by Lj,^ : TqG — > T^G. The smooth subbundle Hi of the tangent bundle TG is said 
horizontal and henceforth denoted by H. We will set V := H2® ■■■ ® Hk and call V the 
vertical subbundle of TG. We shall set v := dimF. One has 

V := codimi/, n = h + v. 

As before, we will assume that dimi/j = hi {i = l,...,k) and that H is generated by 
some basis of left-invariant horizontal vector fields Xh '■= {Xi, This one can be 

completed to a global basis (frame) of left-invariant sections of TG, X. '■= {^i, 
which is graded or adapted to the stratification. This can be done by re-labelling the 
canonical basis {cj : i = 1, ...,n} of g = in a way that it turns out to be adapted to 
the stratification and then by setting 

dx • y 



Xi(x) := L^*e 



Ci [i = L,...,n) 

y=0 



dy 

We shall set ni := hi + ... + hi (no = ho := 0, = n) and Hi = span^{Xi : < i < n/} 
{l = l,...,k). 

Notation 1.8. We shall set Ih := {I,..., hi}, Ih2 := {hi + l,...,n2(= hi + h2)},..., 
Ih^ := {rik^i + l, ...,nk{= n)} , andiv ■= {hi + 1, ...,n} . Moreover, we will use Latin letters 
i,j, k, for indices belonging to Ih and Greek letters a,/3,7, for indices belonging to 
Iv ■ Unless otherwise specified, capital Latin letters I,J,K,..., may denote any generic 
index. We also define the function ord : {l,...,n} — > {l,...,k} by ord(J) := i if, and 
only if, Ui^i < I < Ui {i = l,...,k). 

If p G G and X G we set ^{x,p)it) := expg[tX]{p) (t G M), i.e. '&{x,p) denotes the 
integral curve of X starting from p and it is a 1-parameter sub-group of G. The Lie group 
exponential map is then defined by 

expg : s I — ^ G, expg(X) := expjX](l). 

It turns out that exp^ is an analytic diffeomorphism between q and G whose inverse will 
be denoted by log^. Moreover we have 

dix,p){t)=p*exp^{tX) VtGM. 

From now on we shall fix on G the so-called exponential coordinates of 1st kind, i.e. the 
coordinates associated to the map iogg. 

As for any nilpotent Lie group, the Baker- Campbell-Hausdorff formula (see |16j ) 
uniquely determines the group multiplication • of G, from the "structure" of its own 
Lie algebra q. In fact, one has 

exp,(X) . expg(F) = exp,(X ^ Y) (X, F G s), 

where ^ : g x g — > g is the Baker-Campbell-HausdorfJ product defined by 

X * F = X + r + ^[X, F] + ^[X, [X, Y]] - ^[Y, [X, Y]] + brackets of length > 3. (4) 

Using exponential coordinates, (jlj) implies that the group multiplication • of G is 
polynomial and explicitly computable (see [IS]). Moreover, = expg(0, 0) and the 
inverse of x G G (x = exp (xi, ...,x„)) is x~^ = exp {—xi, — x„). 
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Notation 1.9. Using exponential coordinates for G, every point x = expg{Y2j xjej) G G 
can be regarded as n-tuple x = (xi, ...,Xh,Xh+i, G M". We shall set 

Xh := (xi, Xh) G M^', Xv := (xh+i, a;„) eW {n = h + v). 

Hence x = exp^^Xn ,Xv) = {xh ,Xv)- 

When we endow the horizontal subbundle with a metric Qh = (■, ■)h, we say that G 
has a suh-Riemannian structure. Is important to note that it is always possible to define 
a left-invariant Riemannian metric g = {■,■) in such a way that the frame X. turns out 
to be orthonormal and such that g\H = gn- For this, it is enough to choose a Euclidean 
metric on g = TqG which can be left-translated to the whole tangent bundle and, by this 
way, the direct sum ([3]) becomes an orthogonal direct sum. 

Since for Carnot groups the hypotheses of Chow's Theorem trivially apply, the Carnot- 
Caratheodory distance du associated with gn can be defined as before, and dn makes G 
a complete metric space in which every couple of points can be joined by (at least) one 
-geodesic; see [8], [51]. 

We remind that Carnot groups are homogeneous groups (see [60]), i.e. they have a 1- 
parameter group of automorphisms 5t : G — ^ G (t > 0). Using exponential coordinates, 
we have 5tX = expg{'^j ^^P Xi.Si.) for all x = expg(^^. Xj^.ej^.) G The homogeneous 

dimension of G is the integer Q := Yl^=i ^ coinciding with the Hausdorff dimension of 
(G, (iff) as a metric space (see [H], [M], [35]). 

We remind that the structural constants of the Lie algebra q associated with the (left- 
invariant) frame X. are defined by 

Cfj ■.= {[Xj,Xj],Xr) {I,J,R=l,...,n). 

They satisfy the customary properties: 

• Cfj + Cfj = (skew-symmetry) 

• Ej=i C'j.Ci,, + C'j,,Ci^ + C^jrCI,^ = (Jacobi's identity). 

The stratification hypothesis on the Lie algebra implies the following further property: 

Xi E Hi, Xj G Hh =^ [Xi,Xj] G Hi+rn- (5) 
Therefore, if z G Ihs and j G hr , one has 

C^^O^me Ih^+^ . (6) 
Definition 1.10. Throughout this paper we will use the following notation: 

(t) C% := {Cf^^,,<.l, G MK,y.Km (« ^ In, ); 

(ii) ■■= [C?j]i,j=i,...,n e A<„xn(K) (« e Iv). 
Furthermore, for any Z = Xlag/p. ^aXa E V , we will set 
•^Hcrc, j G {1, ...,k} and ij S Ih^ = {jij-i + 1, ...,nj}. 



12 



M C{Z) := E.,/, z^c-. 

Definition 1.11. Let i G {1,...,A;— 1}. Then the i-th curvature of H is the skew- 
symmetric, bilinear map 

VtH^-H ®H, — > Hi+i, {X ® Y) := [X, Y] mod T 

whenever X E H andY G Hi. By definition ofk-step Carnot group, one has Qh,. (■, ■) = 0. 

Since the bracket map [■, ■] : H ® Hi — > Hi^iii = 1, k — 1) is surjective, the 
definition is well posed. We stress that the 1st curvature Qh{-, ■) := ^h-^ (-, ■) of H is the 
"curvature of a distribution"; see [32], [36], [51] . 

If y G TG, let Y = {Yi, Yfc) be the canonical decomposition of Y with respect to 
the Carnot grading, i.e. Y = Yli=i'P"iO^)y where Vh^ is the orthogonal projection onto 
Hi. Set 

k-l 
i=l 

for X G i/ and Y G TG. Then we have 

Lemma 1.12. Let X & H and Y, Z E TG. Then we have 

(t) {nu{X,Y),Z) = {C„{Z)Y,X)- 

(ii) (^](X,F),Z) = (C(Z)F,X). 
Proof. The proof is an immediate consequence of Definition 11.111 and Definition ll.lOi 

□ 

In the sequel, we will give a definition of connection which recovers the usual definitions 
of Riemannian, partial and non-holonomic connections. Classical notions of connection 
(linear, affine or Riemannian) and related topics can be found in [3H], [SH]- 

Partial connections was defined by Z. Ge in [32]; see also [36] and [H]. Non-holonomic 
connections were used by E. Cartan in his studies on non-holonomic mechanics and by 
the Russian school; see the survey by Vershik and Gershkovich, [62] . 

Definition 1.13. Let N be a C°° smooth manifold and let vr^ : E — > N, : F — > N be 
smooth subbundles of TN . An i?-connection V*-^'^'' on F is a rule which assigns to each 
vector field X G C~(X,E) an R-Unear transformation V^^'^^ ■ C°°(X,F) — > C~(X,F) 
such that 

V X, F G C~(X, E)\/ Z e C°°(X, F) 
W f,ge 

vx, r G c°°(x,E) v/gc°^(x). 

If E = F we shall set := V*-'^''^'' and call an ^'-connection. Any such connection 
will be called a partial connection of TN. If E = TN, then V*-"^^''^-' is called a 



'0 fX+gY^ 



fV^^'^'Z + g\/'y''>Z 



{E,F) , 



U)Vf^^fY =fV'^''>Y + iXf)Y 



{E,F) 
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non-holonomic F-connectioiJ^. If E has a positive definite inner product gs , then an 
E-connection is said metric preserving if 

(m) Zg, (X, Y) = g, (V|X, Y) + g, (X, V'^F) WX,Y,Ze C°°(iV, E). 

The torsion Te associated to the E-connection is defined by 

Te{x, y) := v^r - VyX - Ve [X, Y] VX, Y e C°°(X, E), 

where Ve '■ TN — > E denotes the orthogonal projection onto E. An E-connection is 
torsion free if Tb(X, F) = for every X, F e C~(X, E). We shall say that is the 
Levi-Civita E'-connection on E if it is metric preserving and torsion-free. Note that if 
E = TN , terminology and definitions adopted here are the customary ones. In this case, 
we will denote by V the Levi-Civita connection on TN. 

We stress that the difference between the definitions of partial and non-holonomic 
connection is that the latter allows us to covariantly differentiate along any curve of 
N whereas using the first one only curves that are tangent to the subbundle E can be 
considered. 

Definition 1.14. In the sequel, V will denote the unique left-invariant Levi-Civita 
connection on G associated with the fixed left invariant metric g. Moreover, for any 
X,y e X{H) := C°°{G,H), we shall set V^F := PniVxY). We note that V" is a 
partial connection, also called horizontal connection or //-connection. For notational 
convenience, in the sequel we will denote by the same symbol the non-holonomic connection 
onQ, I.e. V" = V^^^'^). 

Definition 1.15. We define the horizontal curvature Hh of the H -connection 'W" to 
be the trilinear map Rh : H x H x H — > H given by 

Rh (X, Y)Z:= V^^ V^Z - ViV^^Z - [y^^^Z, (7) 

where X,Y, Z ^ X{H). In the sequel, R will denote the Riemannian curvature tensor, 

defined by 

R(X, Y)Z := VyVxZ - VxVyZ - V[yx]Z (X, Y,Z e X(G)). 



Remark 1.16. From Definition \1.14\ using the properties of the structural constants of 
any Levi-Civita connection, we get that the horizontal connection is flat, i.e. 

Vix,=o ii,jeln). 

Note that turns out to be compatible with the sub-Riemannian metric gn , i.e. 

X(F, Z)„ = (V|F, Z)n + {Y, V\Z)n M X,Y,Z e ?i{H). 

This follows immediately from the very definition of 'W" using the compatibility property 
of the left-invariant Levi-Civita connection V with respect to the Riemannian metric g. 
Furthermore, is torsion-free, i.e. 

v|r - VyX - [X, r] = vx,FeX(i/). 

In particular, it turns out that the horizontal curvature R^ is identically zero. 



"'This definition recovers the usual one of "vector bundle connection" (see [l^) where instead of a 
generic vector bundle vr : F — > N we make use of a subbundle of the tangent bundle. 
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Definition 1.17. If ip & C°°(G) we define the horizontal gradient of ip, gradfjip, as 
the unique horizontal vector field such that 

{grad„i), X)h= dilj{X) = VX G X{H). 

Later on, we will denote by Jh the Jacobian matrix of a vector-valued function, computed 
with respect to any given orthonormal frame = {ri, ...,Thi} for H. 

We shall now define the left -invariant co- frame ui := {uj : I = 1, rzl dual to X (with 
respect to to the left invariant metric g). The left-invariant 1-formsU are uniquely 
determined by the condition: 

uJi{Xj) = {Xj, Xj) = 6f (Kroneker) (/, J = 1, n). 

The Cartan's structure equations for the left-invariant co-frame lj are given by: 

n n 

(I) duj = ^uJij Aujj, (II) dujK = ^ujjlAujlk -^JK (/, J, K = l,...,n), 

J=l L=l 

where ujij{X) = (VxXi, Xj) are the connection 1-forms for while VLjk are the curvature 
2-forms, defined by 

Qjk{X, Y) = uk{R{X, Y)Xj) (X, Y e X{G)). 

For what concerns the theory of connections on Lie group and left-invariant differential 
forms, see [38]. Moreover, for many topics about the geometry of nilpotent Lie groups 
equipped with a left-invariant connection, see [17]; for the Carnot case see [501 [51]. 

Remark 1.18. We have 

n 

Vx.Xj = - Y^{Cfj - C'j^ + Cij)XR (/, J = 1, ...,n). (8) 

R=l 

This formula and condition ([5]) are needed to make explicit computations in terms of the 
structural constants. For instance, from ^ it follows that the 1st structure equation for 
the coframe u, becomes 

dujii = ~ J2 ^LJ "^i ^ (9) 

l<7,J<ni_i 

where R E Ih^ = {j : < j < n,} and ^ = 1, k. 

In the sequel we will need the notion of covariant derivative along a path; we refer the 
reader to [11] for a detailed introduction. 

^That is, L'^uJi = loj for every a; G G. 
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Definition 1.19. Let 7 : [a, 6] C M — > G be a path and let X : [a,b] cR — > TG, 
X — Y2i^i{t)-^i{'y)! ^ vector field along 7. Then the covEiriant derivative of X 
along 7, denoted by Vt, is defined as 

n n 

WtX:^V^X = J2{ii+ E rSx(7)07x}^L(7), 

7=1 J,K=l 

where V denotes the Levi-Civita connection on G and 

^JK {^x^Xj, Xi) = u;ji{Xk) (/, J,K^l,...,n) 

are the Christoffel symbols ofV, with respect to the left invariant frame X = {-'^i, Xn\ 
on G. 

We end this section with some important examples. 

Example 1.20 (Heisenberg group H^). Let \)i :— TqIHI^ = M? denote the Lie algebra of 
the Heisenberg group H^, that is the most simple example of2-step Carnot group. Its Lie 
algebra [)i satisfies: 

[ei,ei] = 63 

and all other commutators vanish. We have f)i = i/ © Mea where H = spanj^jei, 62}. In 
particular, the 2nd layer of the grading Res is the center of the Lie algebra These 
conditions determine the group law • via the Baker-Campbell-HausdorfJ formula. Indeed, 
ifx^ expg{Yfi^-^ XiBi), y = expg(^J^^ yiBi) e one has 

X • y = expg{xi + yi, X2 + ^2, 2^3 + ^3 + - X2yi)). 

The standard frame of orthonormal left invariant vector fields for is given, using 
exponential coordinates, by 

X2 

Xi{x) := Lx^ci = Ci - y 63; 

X2{x) := L^*e2 = e2 + y e3; 

Xs{x) := Lx*e3 = e3. 

Example 1.21 (Heisenberg group H"). Let f)„ := TqW = IR^^+i denote the Lie algebra 
of the Heisenberg group H", that is the most important example of2-step Carnot group. 
Its Lie algebra f)„ is defined by the following rules: 

h 

[ei, ei+i] = e2n+i for every i^2k-l, = 1, n = - 

and all other commutators vanish. One has i)n — H (B Me2n+i, where 

H — span][j{ej : i — 1, 2n}. 

The 2nd layer of the grading ]Re2n+i is the center of The above conditions uniquely 
determine the group law • via the Baker- Campbell-Hausdorff formula. More precisely, if 

X = (^Wg{YA=t^ ^i^i): y = ^^PgiTdll^ Vi^i) ^ > then 

1 " 

X»y = expg^Xi + yi, X2n + Vln, X2n+\ + 2/2n+l + ^^{Xlk-\y2k " X2ky2k-l)^ ■ 

k=l 
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If i G {1, ...,n}, the standard frame of orthonormal left invariant vector fields for W is 
given, using exponential coordinates, by 

-^21-1(3;) := Lx^e2i-l = e2i-l ^^2n+l] 

V ( \ — T — I ^2i-l 
-^2i\X) :— ljx*G2i — G2i -\ — G2n+l] 

X2n+l{x) := Lx^e2n+1 = e2n+l- 

Example 1.22 (2-step Carnot groups). We have g = H (B H2, where H2 is the center of 
Q. Furthermore, we have the following rules: 



for every i, j E Ih = {1, hi} 



and all other commutators vanish. Let x = expg(^"^-^ Xjej), y = expg(^"^-^ ?/jej) G G. 
Then 



x»y = expg(x + y-^ ^ {ChX, y)ec 



The standard frame of orthonormal left invariant vector fields for G is given by 
Xi{x) = L^^ei = ei ^ {CHX,ei)ea {i e Ih) 



a&lHo 



■^a{x) Lx*^a 

Remind that C% = [C^] ij^ig ; see Definition \l_AOi 



{a G Ih2 



Example 1.23 (Engel group K^). The Engel group is the simpler example of a 3-step 
Carnot group. Its Lie algebra e is 4- dimensional and is defined by the following rules: 

[61,62] = 63, [61,63] = 64 

and all other commutators vanish. We have c = H 
and the center of the Lie algebra e is M.e4. Therefore Cf, 



3 ©IR64, where H = spanig{6i, 62} 
1 



-1 



and 
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-10 





The group law • is given, for x = expg{Y,t=i ^iei), y = expg(^^^^ j/j6i) G by 



X 



y = expg(^x + y- ^{Clx,y)e3 - (^^{C^x,y) + ,?///) ((^^63, {x-y))^e^. 
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The standard frame of orthonormal left invariant vector fields for ¥} is given by 

Xi{x) := La;*ei = - ^((:7|x//,ei)e3 - (^^(C^Xjei) + ^(C'|xi,,ei)(C''^e3,x)je4 = 1, 

Xj,{x) := Lj,*e3 = 63 - ^(C^x, 63)64 
^4(3;) := ij^*64 = 64. 

Example 1.24 (3-st6p Carnot groups). We have that q = H (B H2 (B H3 where is the 
center of g. In order to describe q, we make use of the matrices of the structure constants 
ofd, i-e. 

Ch = [C°]i,jeiH (aelH^) 
= [C?j]i,j=i,...,n i(3elH,). 

The group law • is given, for x = exp^{YA=i ^i^i), y = expg{YH=i Vi^i) ^ by 
x»y = expJx + y ^ {C'^XH,yH)ea 

The standard frame of orthonormal left invariant vector fields for G is given by 



aelH2 /3e/H3 <^&Ih2 



e/3 



X^{x) = Lx^e^ = ^ {C'^x,e^)ep {'y e IH2); 



W6 r6mind that, for any z = expg{J2i ^i^i) G G W6 usually S6t Zh := (-Zi, Zh) G M.^; 
see Notation 11.91 

1.2 Hypersurfaces: some basic facts 

W6 now r6mind som6 basic facts about Iiyp6rsurfac6s which will be needed in the sequel. 
This material can be found in [501 IM] . 

Later on, and 5™ will denote, respectively, the usual and the spherical Hausdorff 
measures associated with the CC-distance. The (left-invariant) Riemannian volume form 
on G is defined as a'^ := A.'Li^i G A"(rG). 

Remark 1.25. By integrating cr^ we obtain a measure voV^, which is the Haar measure of 
G. Since the determinant of Lx^, is equal to 1, this measure equals the measure induced on 
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G by the push-forward of the n- dimensional Lehesgue measure £" on R" = g. Moreover, 
up to a constant multiple, vol'^ equals the Q-dimensional HausdorfJ measure Ti^^ on G. 
This follows because they are both Haar measures for the group and therefore equal up to 
a constant; see [5^ . Here we assume this constant equal to 1. 

In the study of hypersurfaces of Carnot groups we have to introduce the notion of 
characteristic point. 

Definition 1.26. If S G G is a -smooth (r = 1, oo) hypersurface, we say that S is 
characteristic at x E S if dim Hr^ = dim(i4 fl T^S) or, equivalently, if C T^^S. The 
characteristic set of S is denoted by Cs, i-e. 

Cs:={xeS : dimH^ = dim{H^ f] T^S)}. 

A hypersurface S (Z G, oriented by its unit normal vector u, is non- characteristic if, 
and only if, the horizontal subbundle H is transversal to S {H &\ TS). We have then 

rh T^S ^ Vhv{x) ^ ^ 3X G X(i/) : (X(x), v{x)) ^0 (x G S), 

where Vh '■ TG — > H denotes the orthogonal projection onto H . 

Remark 1.27 (Hausdorff measure of C5; see |45j). If S d G is a C^-smooth hypersurface, 
r > 1 then the Q — 1-dimensional Hausdorff measure associated with dn of Cs is zero, i.e. 

ng'\Cs) = o. 

Remark 1.28 (Riemannian measure on hypersurfaces). Let S G G be a -smooth 
hypersurface and let v denote the unit normal vector along S . By definition, the n — 1- 
dimensional Riemannian measure along S is given by 

crr'L5:=(z/JaS)|5, (10) 

where J denotes the "contraction" , or interior product, of a differential forn^. 

Since we shall study smooth hypersurfaces, instead of the usual weak definition of 
if-perimeter measure (see [3], [ID], [10], [T7], [271 IMl |29] , [31]) we now introduce a (ra — 1)- 
differential form which, by integration, coincides with the i/-perimeter measure. 

Definition 1.29 (crS~^-measure on hypersurfaces). Let S C G be a -smooth non- 
characteristic hypersurface and let us denote by v its unit normal vector. We will call 
i/ -normal along S, the normalized projection onto H ofv, i.e. 

VhV 



We then define the {n — 1) -dimensional measure (t2~ along S to be the measure associated 
with the (n — 1)- differential form a^"^ G A"~^(r5') given by the contraction of the volume 
form of G with the horizontal unit normal v^, i.e. 

aT'^S:={u„Jal)\s. (11) 

If we allow S to have characteristic points, we may trivially extend the definition of a]^''^ 
by setting aS^^L = 0. We stress that a^-^\_S = \Vhu\ ■ a'^'^ L S. 

^The linear map J : K^{T<G,) A'^-i(TG) is defined, for X ^ TG and uj^ G A'=(TG), by 
(XJw'=)(ri,...,n_i) :=c^'=(X,ri,...,rfe_i); see [3H], M- 
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Definition 1.30. If is the horizontal unit normal along S, at each regular point x G 
S \ Cs one has that = {i'jj)x © HxS, where we have set 

H^S '.= Hx n TxS. 

We call HxS the horizontal tangent space at x along S. We define in the obvious way 
the associated subbundles HS{c TS) and v^S, called, respectively, horizontal tangent 
bundle and horizontal normal bundle of S. Moreover we shall set: 

(^) ^ '■= \T^\ = i^H,^), 'Where 



(ill) Ch{zu) := Eae/«, ^aC"^] 
(iv) C{w) := E.6/,^a^" 

2 The CC-distance in Carnot groups 

2.1 On normal and abnormal CC-geodesics in Carnot groups 

Below we shall introduce the main notions about normal and abnormal geodesies in the 
setting of Carnot groups and we shall explicitly write down the associated equations. Our 
approach is the Hamiltonian one and we follow that given in [34J (see also [55]). However, 
at the end of this section, we shall also briefly recall and discuss the Lagrangian point of 
view, which will turn out to be useful in the sequel. 
From now on we shall set 

h := hi = dimH. 

As already recalled in the previous section, the CC-metric dn measures the distance 
between two given points p and q by minimizing the length of all (absolutely continuous) 
horizontal curves (i.e. tangent to the horizontal subbundle H C TG) joining p and 
q. Thus we have to study minimizers for this metric. A minimizer is any absolutely 
continuous horizontal curve 7 : / C M — > G which is such that for every t G / there exists 
e > such that 7 minimizes the length between 7(to) and 7(^1) whenever to?^i belong to 
{t — e,t + e) (Z I . A specific feature of sub-Riemannian geometry is that minimizers can 
be of two different types. However they are not necessarily mutually exclusive. Roughly 
speaking, the minimizers of the first type, called normal, are projections of solutions of 
a Hamiltonian system and, in a sense, they generalize the Riemannian situation since, in 
particular, they are differentiable. Minimizers of the second type are called abnormal or 
singular. Although their existence was originally deduced from the Pontrjagin Maximum 
Principle (see [54J and discussion therein), they can also be defined as projection onto G 
of characteristic curves (in the symplectic sense) of the annihilator of H in the cotangent 
bundle r*G, as we shall see below. 
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For sake of completeness we recall the main definitions of the hamiltonian formalism in 
our particular setting. A Hamiltonian is a function Ti : T*G — > M where the cotangent 
bundle T*G is the phase space. The exponential coordinates x = exp^^xi, ...,a;„), which 
have been fixed on the whole G, induce fiber coordinates on T*G by expanding any 
arbitrary covector p G T*G in terms of the coordinate covector fields {dxi, ...,dxn}, 
i.e. p = Yld=i'Pidxi. The 2n functions {x,p) = (xi, are said canonical 
coordinates on the phase space and the 1-form 6 = X]/=i Pidxj is the tautological one form 
on T*G (which is actually independent of the choice of coordinates on G). The canonical 
symplectic form on r*G is, by definition, the non-degenerate 2-form u = —dQ. Now, 
given a function if, uj uniquely determines a vector field X-^ satisfying dTi = uj^Xt-c, ■), 
which is called the Hamiltonian vector field for Ti, (or symplectic gradient of Ti). The 
Hamilton's equations for a smooth Hamiltonian Ti are the O.D.E.'s for the integral curves 
of X-H. In canonical coordinates, they are given by 

XI = Pi = -jr- a = i,->^)- (12) 

opi OXi 
We remind that the momentum function Py '■ T*G — > M is defined by 

Py{x,p) ■■= piY{x)). 

In the sequel pi will denote the momentum function associated to the J-th coordinate 
vector fielc|llc^/c^a;/. Hence, in canonical coordinates one has Py{x,p) = Yll=iPi^ii^)^ 
where Y = Yl^=i Yid/dxi. 

We now start with the derivation of the CC-geodesic equations by defining the suh- 
Riemannian Hamiltonian (or kinetic energy) Tisiz ■ To this aim let us set 

Pi := Px, for every i e Ih = {1, h}, 

to denote the momentum functions associated with a orthonormal (left-invariant) moving 
frame = {Xi, ...,Xh} for H and note that, if Xi{x) = ^"=i(Xj(x))/e/, one has 

n 

Pi{x,p) = PxA^^P) = '^iXii^))iPi for every ie Ih. 



1=1 



Notation 2.1. We shall set 

• Ph := Y.i&lH Pi^i' 

• Pv ■■= Y,a€lv -^a^a, 

^We are using the notation: 



d/dxi=ci^{0,..., ^ ,...,0) (/ = l,...,7l). 

I — th place 
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Definition 2.2. The sub-Riemannian Hamiltonian is defined by 

The Hamiltonian equations 0121) associated with the sub-Riemannian Hamiltonian Tisiz 
are called the normal geodesies equations. A normal curve is the projection onto G 
of a solution of the normal geodesies equations. 

The following result is well-known and a proof can be found in [54]. 

Theorem 2.3. Every sufficiently short arc of a normal curve x G G is a minimizing 
CC-geodesic. Moreover x is the unique minimizing CC-geodesic connecting its endpoints. 

The normal CC-geodesic equations are given by the following system: 



(13) 

see Definition 11.101 for the notation C". Note that the first equation express the fact that 
the normal curve is horizontal. The second equation can be deduced by noting that, for 
any function / : T*G — > M and for any solution x : / C M — > G of the Hamiltonian 
equations, one has 

|/(7W) = {/,^«.}, 

where {■, ■} denotes Poisson bracket. In particular P/ = {Pi,Hsr} for every / = 1, ...,n. 
Therefore, it follows that 

Pi = {Pi, \ E ^'(^'P)} = \ Y^^Pi^P^ = " E E 

jG/h ialn *S/_f/ QS-fi/ 

The last computation yields ([T3l) . 

Below we shall discuss some other features of the normal geodesic equations but let 
us first introduce the following (see |3l], [51]): 

Definition 2.4. An abnormal curve is a horizontal curve which is the projection onto 
G of an absolutely continuous curve in the annihilator H-^ C T*G of H , with square 
integrable derivative, which does not intersect the zero section of H-^ and whose derivative, 
whenever it exists, is in the kernel of the canonical symplectic form restricted to H^ . An 
abnormal minimizer is an abnormal curve which is a minimizer. A strictly abnormal 
curve -resp. minimizer- is an abnormal curve -resp. minimizer- which is not normal. 

For a deep study on abnormal extremals and related problems in sub-Riemannian 
geometry, see pQ, [2] and references therein. 

The equations for abnormal curves in our Carnot setting can explicitly be derived 
along the lines of |34i (see also [55l IS]). More precisely, they are given by 

I^ae/v PoChXh = 

(14) 



Pv - - Yla&Iv PaC^'X 



Xv = 
Ph =0. 
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Unlike the system f|T3|) . the equations of the system (fT6!l are mixed algebraic-differential 
equations and they cannot be expressed as O.D.E.'s. Note also that abnormal curves only 
depends on H and not on the metric. 

We remind that, for any x = expg{J2i ^i^i) G G we set Xh := {xi, ...,Xh) E = H 
and Xv ■= {xh+i, ■■■,Xn) eW" = V; see Notation 11.91 

Definition 2.5. According to Definition M.lCK we shall set: 
(i) Ch{Pv) := Yl,a&iv ^^Ch; 
(n) C{P.):= EaeiyP^C""- 
Using this notation, (|T3l) and (|T6l) can be rewritten more compactly as follows: 



X = Ph 

(Normal Equations) • n(u w ^^^^ 

P = — C \Pv )Ph ] 



(Abnormal Equations) 



Ch{Pv)Xh = 

Py = -C{Pv)X 

Xv = 
Ph = 0. 



(16) 



In the next examples we shall write down the normal CC-geodesic equations. Some of 
them will be studied in greater detail in the sequel. 



Example 2.6 (Heisenberg group H^). We get 

X = Ph 

Ph = — C H {^Pz)Ph 
P3 = 0. 



More explicitly, one has Ph = —P3 



1 
-1 



Ph , where Ph 



Pi 
P2 



Example 2.7 (Heisenberg group IP). We get 

X = Ph 

Ph = —Ch {P2n+l)PH 
P2n+1 = 0. 



We have 



Ph = -P, 



2n+l 



10 
-10 
1 
0-10 








1 
-1 



Ph., 



Ph 



Pi 
P2 
P3 
Pa 

P2n-l 
P2n 



(17) 
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Example 2.8 (2-step Carnot groups). We get 



X = Ph 

Ph = -C„{Ph,)Ph (19) 
Ph, = 0, 



where Cn{PH,) = Y.o.^,^^Pc.C'^^- 
Example 2.9 (Engel group ¥}). We have 



X = Ph 

P, = -{P,Cl+P,C')Pn 

Pg = -P^C^P„ ^^^^ 
Pi = 0. 

Example 2.10 (3-step Carnot groups). For the general 3-step case from ( |T5l) we get 



(21) 



X — Ph 

Ph = -Ch{Ph,)Ph 
Ph, = -C{Ph,)P„ 

A3 = 0, 

where Ph, := Y^aein, -^a^"' -^^3 := Y^aein^ ^^^^ remind that 

Ch{Ph,) = J2 Po.C^^ C{Ph,) = J2 Pc.C''- 

2.2 CC-normal geodesies, variational formulae and Jacobi fields 

In this section we shall again discuss normal CC-geodesics and their minimizing properties. 
We first recall some results which can be found in [62] . 

If we look for minimizing geodesies joining two (fixed) points x, y E G, we are solving 
the minimization problem (with fixed endpoints): 

f' 

min / \xH\dt, x(t) G H^^t) Vt G [a,b], x = x{a), y = x{h). 

J a 

Notation 2.11. If Z E X(G), then Zh-, Zv denote the orthogonal projection of Z onto 
H and V , respectively. 

The Euler-Lagrange equations for this problem, written with respect to the global left 
invariant frame X = {Xi, ar^ 

- Lj; = ^ {Pa^a - PaixH J duJa))^ 



^Hereafter the symbols b and will be used to denote the so-called musical isomorphisms between 
vectors and co- vectors, defined in an obvious way with the help of the metric; see [H]. For the definition 
of the "contraction operator" _l , see footnote [HI 
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where L{t,x,x) = \xh \ and Pa {a G Iv) denote the a-th Lagrange multiplier, see [62] . 
The horizontahty for the curve x is expressed by the equations co'q(x) = (a G Iv). 
These equations can be rewritten in invariant form by replacing the Lagrangian L of the 
unconditional problem by Lsjz := L + ^^g/^ PaUJa{x). Using the Levi-Civita connection 
related to the fixed (left invariant) Riemannian metric on G, the normal CC-geodesic 
equations are given (see ^62j) by 

VfXH + Y^aelv {Pa^a - PaiXn J duJa)) = ^22) 

uja{x) = (a G Ji/). 

Note that the second equation is again the horizontahty condition which can equivalently 
be expressed by Xv = 0. Setting Pv := J2aGiv C'(Pv) := '^^gIv 

system f l37|) can be rewritten as 

+ Pv +C{Pv)xh =0 
Xv = 0. ^^^^ 

To see this it is enough to express the right-hand side of the first equation in ( 1371) using 
formula ([9]) for the exterior derivative dua (a G Iv)- Note that the Lagrangian multiplier 
Pv can be regarded as a curve in the vertical subbundle V. Obviously, (ITS!) and (1231) 
coincide. This immediately follows by explicitly calculating the first equation in (123|) . 
More precisely, as in the Riemannian cas^, it turns out that 



L=l i,j&lH 



where F^- := {V x^Xj, Xl) denote Christoffel Symbols. Using formula ([8]) together with 
condition on the Carnot structure constants, we get that F^- = |Cj^ = whenever 



L E Ih {i,j & Ih). Moreover, by skew-symmetry of the structure constants C[<, for everj'^'^' 



L E Iv we get that J^ijeijj ^ij^i^j — 0- Therefore Xh = VtXu- Setting Ph := Xh and 
projecting (l23ll onto H and F, respectively, one gets 



X — Xh — J^H 

Xi/ = 

Ph =Xh = VtXH = -Ch{Pv)Ph 
Pv = -C{Pv)Ph 



(24) 



which is the projected form of (1151) . as we wished to prove. 

As already said in the previous discussion, the normal CC-geodesic equations can be 
directly deduced by minimizing the constrained Lagrangian 



I/snyti X, xj — \Xh \ + {Pv,x) = \X„\+ Pa 



^UJa{X). 



^See |TT], [35], [32] for a classical setting, or for a discussion of geodesies equations for Nonholo- 
nomic geometries. 



Actually, by using again ©, C^^ can be different from only if L belong to Ii 



«2 ■ 
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This is equivalent to compute the so-called 1st variation of the functional 



Isiz(x) :-- 



Lsiz {t,x,x)dt. 



For this reason and in order to explicitly computing the 2nd variation of Isiz{x), we need 
preliminarily the next: 

Definition 2.12. If x : [a, b] — > G is a smooth path, we call variation of x any smooth 
mapping 'd : [a,6]x] — eo,eo[ — ^ G , eo > 0, for which x{t) = 0) for all t G [a,b]. 
We say that ■(9 fixes endpoints if x{a) = ^9(0,5), x{b) = 'd{b,s) for all s e] — eo,eo[- 
In this case we say that d is a homotopy of x. Moreover, if for every s g] — eo,eo[ 
the path 'ds : [a, 6] — > G, 'dsif) = i^(t,s), is a normal CC-geodesic, we say that -d is a 
sub-Riemannian or CC-geodesic variation of x. 

In the sequel we shall write dt,ds for ■d^{dt),^^{ds) respectively, and denote covari- 
ant differentiation of vector fields along with respect to dt, ds by Vt, respectively. 
Furthermore, let us recall the following general identities: 



where R is the Riemannian curvature tensor of G. 



(25) 
(26) 



Proposition 2.13 (1st variation of Isr.{x)). Let us assume that x : [a,b] — > G is a 
differentiable path and let ^ : [a,6]x] — eo,eo[ — ^ G be a differentiable variation of x. 
Then 



S^Isii 

J dtx 
^\dtx\ 



d 
ds 



s=0 



+ dtPv + C{Pv)dtX 



\dtx 

{dsPv\s=o,dtx) \dt. 



(27) 



Let us assume that \dtx\ = 1 and let ^ be a homotopy of x. Setting 

Y{t) := ds^{t, 0), Qv{t) := dsPv {t, 0) 

we get that 

S^Isjz = 



{Qv,dtx) 



Vtdtx + dtPv + C{Pv)dtxj^yt. 



{2i 



Note that (123!) immediately follows from the previous result once we require that the 
path X is an extremal of the functional Isr.{x). Indeed, in such a case, the fist term under 
integral sign in (|28l) must vanish for every F, which gives the first equation in (l23ll . The 
same procedure implies that the second term under integral sign vanishes for every Qv, 
which in turn implies the horizontality condition for x. 
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Proof of Proposition \2.13\ . First we have to compute 

ds [ {\dt^H\ + {Pv,dt^)\dt = ds [ \dt^H\dt + dJ {Pv,dt'&)dt. 

J a J a J a 



= :Ao 



Actually, the first term Ai can be directly deduced from the Riemannian computation 
(see, for instance, |TT], Theorem 2.3) and, with our notation, we get that 



So we have to compute the other term. We have 
A2 = ds [ {Pv,dt^)dt= I ds{Pv,dt^)dt 

J a J a 

= ^ {{VsPv,dt^) + {Pv,\/sm]dt 

= {dsPv + Vds&Pv , dt^) + {Pv , Vtds^)yt 

= [ {{dsPv,dt^) + {V9,^Pv,dt'd) + dt{Pv,ds'&) - {VtPv,ds^)]dt 

J a 

= ^ {{dsPv,dt^) + {V9,^Pv,dt^)+dt{Pv,ds^)-{dtPv+V9,^Pv,ds^)]dt 

= {dJ,Pv)\l + j { {dsPy , dt^) - {dtPv , ds'd) + {Va^^Pv , m - {VtPv , ds^l ]dt. 

" =:B 

We claim that B = {[ds'&, dt'&], Pv)- To prove this claim we may proceed by using a well- 
known formula for the Riemannian connection which can be found in [11] (see formula 
1.29, p. 15). This way we get 

B = {V9,^Pv,dt^)-{VtPv,ds^) 

= ^[{[ds^,Pv],dt^) - {[Pv,dt^],ds^) + {[dt^,ds^]] 

- ^^^{[dt^,Py],ds^) ~ {[Pv,ds^],dt^) + { [ds^, dt^],Pv)] 

= {[ds^,dt^],Pv). 

Moreover, it is easy to show that {[dg^, dt'd], Pv) = —{C{Pv)dtX, dg^). Indeed one has 

n n 

{[ds^,dt^],Pv)= $^(5.^9)/(ai^9)jP„([X,,Xj],X,)= 5^ C;^(9.^)/(5t^)jPa 

I,J=la&Iv I,J=laeIv 

and the claim immediately follows by using (ii) of Definition 12.51 So we have shown that 
A2 = {ds^,Pv)\l + {{dsPv,dt'&) - {{dtPv+CiPv)dt^),ds^)yt. (29) 
By adding Ai and A2 (1271) claim easily follows. Finally, (l28ll follows from (p7|) . □ 
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Theorem 2.14 (2nd Variation of Isr.{x)). Under the notation of Proposition [KTB let x 
be a normal CC- geodesic satisfying \dtx\ = 1, and assume that -d is a homotopy of x. 
Moreover, set Y(t) := ds^{t,0) and Qvif) := dsPv{t,Q). Then we have 



5Hs 



Vf + C{Pv ) {VtY„ + [r, dtx\) + [F, dtx] + R(9tx, Y)dtx\ ) } 



>dt 



(30) 



Proof. Since, by hypothesis, x is a normal CC-geodesic and is a homotopy of x, we may 
start by considering the following identity: 



d 
ds 



Isizi^s) = - ^ {{ds^, [Vtdt^H + dtPv + CiPv)dt^Hj ) - {dsPv,dt^)yt. (31) 



This identity can easily be deduced by the proof of Proposition I2.13[ So we have 



ds 



Vtdt^H + dtPv + C{Pv)dt^H] ) - {dsPv,dt^)yt 



a " s. 



Vtdt^H +dtPv +C{Pv)dt^. 
Vtdt^H +dtPv +C{Pv)dtd„ 



ds{dsPv,dt^)\dt 



Vtdt^H + dtPv + C{Pv)dtdH\ ) - ds{dsPv,dt^)] dt. 



■.A2 



:A3 



Now it is obvious that, at s = 0, one has = because x is assumed to be a normal 
CC-geodesic. 



Step 1. (Computation of A2) We have 
Vs{Vtdt^„ + dtPv + C{Pv)dt^„) 



s=o = "^I'^Yn + [Y,dtx] + R{dtX,Y)dtX + VtQv 
+CiQv)dtx + C{Pv){VtY„ + [Y,dtx]). 

Proof. This computation generalizes the classical deduction of Jacobi's equation. We have 

B := Vsi^tdt^H + dtPv + C{Pv)dt^H) = VsVtdtdH +VsdtPv + Vs{C{Pv)dt'dH) . 



■Bi 



■■B2 



■■B3 



The first term Bi can be computed in analogy with the Riemannian case; see [20], p. 111. 
More precisely, by fl2Bl) we get 



B, = VsVt{dt^)H = VtVs{dt^)H - K{d,§,dt€){dtd)H. 
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At this point we have to compute the term Vs{dt'&)H- Note that this must be done in a 
different way with respect to the Riemannian case because the term dt'& is not horizontal, 
a priori, and so it may be dt^ 7^ dt'&H ■ We claim that 

Vsidt^)H =Vtids^)H + [ds^,idt^)H]. (32) 

Indeed we have 

Vs{dt^)H = Vs{V„{dtd))=ds{V„{dtd)) + Vo^^V„{dtd) 

n 
n 

a&Iv 

and fl32|) follows. Therefore we obtain 

Bi = VsVtidt^)H = VtVtids^)H + [ds^, + R{dt^,ds^){dt^)H. (33) 

Remark 2.15. In what above, we have used the identity 

ds{V„{dt'd))=dt{VnidJ)). 

In order to prove this identity, we may argue componentwise. More precisely, for i Ih 
we compute 

dt{dA Xi) = {m^, X,) + {ds^, dtX,) = {dtdj, X,) + {dj, JXA^). 

By subtracting the second identity from the first one, we easily get that 

Ai ■= ds{dt^,X.^ - dt{ds^,Xi) = {dt'd, [jXi - {JX>)''\d,d)). 
From the last relation we infer that 

n 

Ai=Y ([^^i - {JX,f']Xi,Xj){ds^)i{dt'd)j. 
i,j=i 

Now since 

{JX,Xj,Xj) = -{VxrXj,X,), 

{{JX,)^'Xj,Xj) = {{JX.,f'Xj,Xi) = ~{VxjXj,X,), 
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we get that 

n n 

/or every i & Ih , which implies the claim. 

Remark 2.16. // the variation d is a CC-geodesic variation of x, the term Bi can be 
computed as in the Riemannian case; see fWf . p. 111. Indeed, in such a case we have 
dt'd = Vn^dt'd) for every s g] — eo,eo[. Using fl25l) and fl26l) yields 

= VtVt{ds^) + R{dt^,d,m^- (34) 
This can also be deduced from fl33l) . by noting that 

[d.A {dt^)H] = [ds^,dt^] = A[ds,dt] = 0, 

because [ds, dt] =0. 

The term B2 can be computed by means of fl25|) and we immediately get that 

B2 = VsdtPv =VtdsPv. (35) 
Analogously, the term B^, can be computed by means of fl2B]) and (152]) . as follows: 

= C{dsPv)dt^H + C{Pv)VA^H = C{dsPv)dtd„ + C{Pv){Vt{ds^)H + {dtd)„]). 

Finally, by adding the terms Bi, B2 and -B3, we get 

B = VtVt{ds^)H + [ds^, imn] + RidA ds^)idt^)H + VtdsPv 
+C{dsPv)dt^H + C{Pv){Vtm)„ + [ds^, {dt^)H]). 

The thesis follows by substituting Y{t) = ds'd(t,0) and Qv(t) = dsPv(t,0) into the last 
expression. □ 

Step 2. We have { jj^ A^dt} \^^^ = - J^^ (F, [dtQv + C{Qv)dtx] )dt. 

Proof. By arguing as in Proposition 12. 13[ we get that 

j' Asdt= {dsPv,ds^)\l + J' {(^^f,dt^) - (9,^, [Cid,Pv)dtX + mPv])]dt. 

Now we consider this expression at s = and we make the substitutions: dt'&\s=o = dtx, 
ds^\s=o = Y and dgPv \s=o = Qv ■ By hypothesis, is a homotopy. Hence 



{Qv,Y)\^ = {dsPv\s=o,ds^\s=o)\, = 0- 
Set s = 0. Then (^g^, = {^-^^\ ^^q, dtx) = and the thesis follows. □ 
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Using Step 1 and Step 2, we get that 



Y, 



5Hs 



S^Isiz = = - ^ {(^' (V^^H + [Y,dtx] + R{dtX,Y)dtX 

+VtQv + CiQv)dtx + CiPv){VtY„ + [Y,dtx]) + {dtQv + CiQv)dtx)^^yt. 
Since VtQv = dtQv + VdtxQv , it follows that 

+2(vtgv- +c(gv)9tx- ^Va^.Q;^) +c(Pv)(VtrH + [Y,dtx\)\)yt. 

By reordering a bit the last expression we easily get that 

Vf V„ + C{Pv){VtY„ + [r, S^x]) + [F, + R(9tx, 
VtQv +C{Qv)dtX-]^Va,xQv\)]dt. 
Step 3. r/ie following identities hold true: 

{%) {V9,,Qv,Y) = '^{{C{Yy)Qy,dtX) + {C{Qy)dtx,Y))- 

(ii) {C{Qv)dtx,Y) = -{[dtx,YlQv). 
Proof. To prove the first identity we compute 

n 

1 " 



+2 F, 



(36) 



=Oby (6) 



i((C(rv^)Q.,9ix) + (C(Q.)5iX,F) 



which proves (i), while (ii) follows from Definition 12.51 and Lemma [1.121 
Using Step 3 and integrating by parts, yields 

r-6 



□ 



VtQv +C{Qv)dtx--V9,.Q^ 



dt 



212^^/ { - {"^tY.Qv) - {[dtx.YlQv) -\[{C{Yv)Qv,dtx) + {C{Qv)dtx,Y) 



>dt 
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^ I - {VtY.Qv) - -^{[dtX.YlQy) - ^-{C{Yv)Qy,dtx)]dt 



Q^ 



VtY + ^[dtX,Y]-^CiYv)dtX 



dt. 



Finally, by using the last expression and (15^ . we get that 



6'h 



Y, 



Vf V„ + C{Pv){VtYH + [Y, dtx]) + [Y, dtx] + R{dtx, Y)dtX 



2{Qv, \s/^Y -'^[Y,dtx\~hj{Yv)dtX^)yt, 



which is equivalent to the thesis. 



□ 



Accordingly with the Hamiltonian theory already discussed in Section 12.11 we set 
Ph '■= dfX. Starting from the second variation formula (I3UI) . it is natural to consider the 
following system of O.D.E.'s: 



Vf V« + C{Pv ) {VtYn + [Y, Ph]) + [Y, Ph] + R{Ph , Y)Ph = 
Vv {VtY -I[Y,Ph]- \C{Yy)PH) = 0. 



(37) 



Corollary 2.17 (2nd derivative of Isr.{x) through CC-geodesic variations). Under the 
notation of Proposition IKT3[ let x be a normal CC-geodesic satisfying \dtx\ = 1, and let 
us assume that: 

(i) d is a homotopy of x; 

(ii) d is a CC-geodesic variation of x. 

Moreover, set Y{t) := ds'd{t,Q) and Qv{t) := 9,Py(t,0). Then 

6^Isiz= j {(g,.,(ViF+[9iX,F]))-(r, (vfV + R(aia:,F)9iX + C(Pv)ViF))}rft. 



(3^ 



Proof. The proof follows the same lines of the general case. Since, by hypothesis, a; is a 
normal CC-geodesic and the CC-geodesic variation d is also a homotopy of x, we may 
start by considering the identity 



d 
ds 



Vtdt^ + dtPv +c{Pv)dtd 



dt. 



(39) 
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Indeed note that dt'd = dt'dn for every (t, s) G [a, h] x] — cq, eo[- This identity can easily be 
deduced, as in the proof of Proposition 12. 13[ by using the hypotheses on x and 'd. Hence 



ds 



a - ^ 



Vtdt^ + dtPv +c{Pv)dt^ 



= :Ai 



Vtdt^ + dtPv +C{Pv)dt^ 



dt. 



Now it is obvious that, at s = 0, one has Ai = 0, because a; is a normal CC-geodesic. So 
we need to prove the following fact: 
Step 1. (Computation of A2) One has 



Vs{Vtdt^ + dtPv+c{Pv)dt^) 



s=0 



= vfV + R{dtx,Y)dtx + VtQv + C{Qv)dtx + C{Pv)VtY. 
Proof. The proof mimics that of Theorem I2.14[ We have 



B := Vs{Vtdtd + dtPv + C{Pv)dtd) = YsVtd£ + YsdtPv +Vs{C{Pv)dt^) . 

■B2 



The first term Bi can be computed exactly as in the Riemannian case (see [20], P-lll or 
the previous Remark [2 .lOp . More precisely, by using (l25l) and (l26l) . we get 



Also B2 can be computed exactly as in the previous proof of Theorem 12.141 (see Step 1) 
and we get that 

B2 = VsdtPv =VtdsPv. 
Analogously, for the term B-j, we get 

i?3 = Ws{C{Pv)dt^)=Vs[ Yl Pc^C'^^t^) = E {dsPaC''dt^ + PaC''Vsdt^) 

= C{dsPv)dt^ + C{Pv)Vsdt^ = C{dsPv)dtd + C{Pv)Vt{dsd). 

Finally, by adding the terms Bi , B2 and -B3 we have 

B = Wt^t{ds^) + R{dt^, ds^){dt^) + WtdsPv + C{dsPv)dt^ + C{Pv)Wt{ds^). 

The thesis follows by substituting Y{t) = dsi!^(t,0) and Qv(t) = dsPv{t,0) into the last 
expression. □ 



33 



By what proved in Step 1 we get that 



S=0 



/b 
(vfV + R{dtx,Y)dtx + VtQv + C{Qv)dtx + C(Pv)VtF))}dt. 

Now, by arguing again as in proof of Theorem 12.141 we obtain 

/ {Y,{VtQv +CiQv)dtx)dt= {Y,Qv)\i- {{WtY+[dtx,Y]),Qv)dt 

J a J a 

and since i? is a homotopy of x we have {Y,Qv)\i = 0. Putting all together, (EHD follows. 

□ 

As before (see (137|) ). by setting Ph := dtx and using (1381) . we deduce the following O.D.E.'s 
system: 

Vf V + R(Ph , Y)Ph + C{Pv)VtY = 

(40) 

Vv{VtY-[Y,PH]) = 0. 

Definition 2.18 (Jacobi equations for normal CC-geodesics). We say that (HOl) represents 
the system of Jacobi equations for normal CC-geodesics. Let x : [0, a] — > G be a normal 
CC-geodesic satisfying \dtx\ = 1. Then we say that a vector field J G X{G) is a Jacobi 
field along the normal CC-geodesic x if and only if J satisfies fHOl) for all t G [0, a]. 

Remark 2.19. An interesting corollary of the second variation formula can be formulated 
by considering only some particular CC-geodesic variations. Indeed, we could consider 
variations of a given normal CC-geodesic x, or, more precisely {x,Pv), through normal 
CC-geodesics such that, for every s g] — eo,eo[ one has Qv = 0. Roughly speaking, the 
variation 'd is chosen in such a way that -ds satisfies fl23l) for each s g] — ^o,^o[, with the 
same Pv • By making this substitution into fl55]) we obtain the following formula: 

I'b 

6^IsR = - / {Y, (vfV + R{dtX,Y)dtX + C{Pv)VtY))dt. 

J a 

We therefore obtain the following 2nd order linear system of O.D.E. 's: 

Vf V + R{Pn , Y)Ph + C{Pv )VtY = 0. (41) 

Definition 2.20 (Jacobi equations with constant Lagrangian multiplier Pv)- We say that 
(j4T|) represent the system of Jacobi equations for normal CC-geodesics having the same 
Lagrangian multiplier Pv . Furthermore, let x : [0, a] — > G be a normal CC-geodesic 
satisfying \dtx\ = 1 and having Lagrangian multiplier Pv . Then we say that a vector field 
J G X{G) is a Jacobi field along the normal CC-geodesic {x, Pv) if, and only if, J satisfies 
dUD for allte [0,a]. 
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As in the classical setting, a Jacobi field is uniquely determined by its initial conditions: 
J(0), VtJ(O). To see this, it is sufficient to develop either (HOl) or (HT!) along the left 
invariant frame X. = {^i, ■■■,Xn} by using the very definition of Vj and the hnearity of 
the curvature tensor R{x, ■)x (remind that X — Xh = Ph{x)). Indeed, one easily deduce 
that (HTj) is a linear system of O.D.E.'s of the 2nd order. Hence, for given initial conditions 
J(0), J(0), there exists a C°° solution of the system defined on [0, a]. Thus there exist 
2n linearly independent Jacobi fields along (a;, Pv). 

Remark 2.21. As in Riemannian Geometry, Jacobi equations for normal CC-geodesics 
and Jacobi equations with constant Lagrangian multiplier Pv -together with the related 
notions of Jacobi fields- are important tools and they can be used to analyze the sub- 
Riemannian exponential map and to perform a precise study of the conjugate and cut loci 
of a point. Nevertheless, we will not pursue this task here. 



2.3 Sub- Riemannian exponential map and CC-spheres 

The main references for this section are [58], [59], [62] . 

Starting from the system (ITSl) for normal CC-geodesics, from the standard O.D.E.'s 
theory we easily obtain existence, uniqueness, regularity and smooth dependence on the 
initial data for small times. So let a;o G G be a fixed point. Then we shall denote by 

exp<^(xo,Po) : [0,r] C R — ^ G, exp^^ (xq, Po)(i) 

the (unique) normal CC-geodesic starting from xq with initial condition P(0) = Pq for 
some fixed vector Pq = Ph{0) + Pv{0) G i^o © Vxq- Here r < rmax(-Po) where rmax(-Po) 
is the maximal time of existence and uniqueness of the solution x(t) = exp^ (xq, -Po)(^) 
of (1151) . Actually, by standard "extendibility" results for O.D.E.'s, it can be shown that 
each solution of system f[T^ is globally defined^ on G, i.e. rmax(-Po) = +oo. 

In order to define a sub- Riemannian equivalent to the ordinary exponential map we 
preliminarily state the following: 

Lemma 2.23 (Homogeneity). If the normal CC-geodesic x{t) = exp <^{xq, Po){t) is de- 
fined on the interval [0,r[, then the normal CC-geodesic exp s,z{xo,aPo){t) , a > 0, is 
defined on the interval [0, ^[ and it turns out that 



exp <^ (xo, aPo)it) = exps^ixo, Po)(at) for every t G 



r 

0,- 

a 



^^More generally, in the setting of sub- Riemannian manifolds, it can be shown that each solution of the 
system of normal CC-geodesic can be continued as long as x{t) remains in the base manifold, so blow-up 
in the dual variable P never occurs. This claim can easily be proved using the following: 

Lemma 2.22 ( 61j). Let x{t) be any normal CC-geodesic defined for t G [0,r[ and assume that x{t) 
remains inside a compact subset of the base manifold. Then x(t) can be extended beyond t = r. 

Proof. The proof can be found in '61], Lemma 4.1. □ 

This result can be used to canonically define a sub- Riemannian exponential map; see [61j . 
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r 

0,- 

a 



Proof. Let {y, Q) : [0, ^[ — > TG be the curve given by {y{t), Q{t)) = {x{at), aP{at)). Now 
we claim that {y, Q) satisfies (fT5|) with ?/(0) = xq and Qo = Q(0) = aPo- Indeed, by the 
very definition of we get that y{t) = ax{at) = aPnlat) and that 

Q{t) = a^P{at) = ~a'C{Pv{at))P„{at) = -C{Qv{t))Q„{t) 

for t G [0, ^[ and the claim follows. So by uniqueness we get, in particular, that 

y{t) = exps,^{xo,aPo){t) = x{at) = exp^^ (xq, Po)(at) for every t G 

and the thesis follows. □ 
Definition 2.24 (Sub-Riemannian exponential map). From now on we shall set 

expsR.{xo,-) ■ TxqG — > G, exp5^(xo,Po) = exp^^ (xq, Po)(l). 

The map expgj^{xQ, ■) is called the sub-Riemannian exponential map at xq G G. 

The sub-Riemannian exponential map parameterizes normal CC-geodesics. Note that 
every minimizing curve connecting xq to a point of G\exp (xq, T^^G) is necessarily strictly 
abnormal. 

The map expg^{xQ, ■) plays in sub-Riemannian geometry a similar role with respect 
to the ordinary exponential map in Riemannian geometry. Nevertheless, there are many 
differences and its structure is much more complicated. An important difference is that 
expsjz{xo, ■) is not a diffeomorphism on any neighborhood of the origin in © Vx^- To 
see this it is enough to choose Pq G Vx^ (i.e. Pff(O) = 0); in this case exp5^(xo, Po) = xq. 
Furthermore, in any arbitrarily small neighborhood of the origin in i/^,, © Vx^ there are 
points Po = Pff(O) -|- Pv(0) with Ph{0) 7^ at which the rank of dexp^^ (xq, Po) is not 
maximal. We shall discuss some of these facts later on. 

We begin by stating the notion of sub-Riemannian wave front. 

Let UH denote the set of all unit horizontal vector of H, i.e. UH <z H = E>^~^ and 
define the map 

H\{0}(B V ^ UH (B V P = T^- 

\Ph I 

The sub-Riemannian wave front Wsiz (xo, r) centered at xq and with radius r > is defined 
as the set of points x(r) = exp^;^ (xo, rPo) = exp^ (xo, Po){r) for Pq G Hx^ © Vxq- In other 
words 

WsR{xo,r) = expsj^{xo,r{UHxo® 14o)). 

Note that in the Riemannian setting the wave front simply coincides with the r-sphere. 
In our case only the inclusion (xo,r) ^ WsR{xo,r) holds true. Actually the structure 
of WsR (xq, t) is very complicated and, in general, the sub-Riemannian wave fronts are not 
manifolds. 

Definition 2.25 (Conjugate locus). A point x G exp5^(xo, Tx^G) is said conjugate to 

Xq if and only if it is a critical value of exp^^ (xo, ■), i.e. there exists Pq G Hxq such that 
X = exp sji{xo, Pq) and dexp sr.{xo, Pq) is not onto. The conjugate locus Conj(xo) of xq 
is then the set of all points conjugate to xq. 
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Notice that, from what we have said above, it turns out that xq G Conj(xo). Moreover, 
by Sard Theorem apphed to exp^(xo,-) one gets that Conj(xo) is a set of (Lebesgue) 
measure zero in G. 

Definition 2.26 (Cut locus). Let us fix xq & G and let us choose a normal CC-geodesic 
x(t) = exp sjz{xo, Po)it) . If t > is sufficiently small, rf// (x(0), = t, i.e. x{[0,t]) is 
a minimizing normal CC-geodesic. Moreover, if for some ti x{[0,ti]) is not minimizing 
the same is true for all t > ti. By continuity, the set of numbers t > such that 
(in (x(0), = t is of the form [0,to] or [0, +oo[. In the first case we say that a;(to) 
is the cut point of xq along x(t) while, in the second case, we say that the cut point 
does not exists. The set Cut(xo) defined as the union of the cut points of xq along all the 
normal CC-geodesics starting from xq is called the cut locus ofxQ. 

2.4 A sub-Riemannian version of the Gauss Lemma 

Remark 2.27 (Working hypothesis). Let G be a k-step Carnot group. Throughout this 
section we shall assume that there are not strictly abnormal minimizers in G; see 
Definition \2.4 

Here below we shall perform an explicit and general computation that will be an 
important tool for the rest of this paper. It is somehow based on the validity of the 
eikonal equation, first proven in |53j . 

We may start by the O.D.E.'s system (fTSll . We assume the solution is parameterized 
by arc-length. Furthermore, let us fix the initial conditions: x(0) = xq, -P(O) = Pjll and 
set dn (x) := dn {x, xq). We also assume the solution to be with unit speed (i.e. \Ph\ = 1). 
All together, these assumptions uniquely determine a normal unit-speed CC-geodesic. In 
particular, we get 

dn {x(t)) = t 

for every t G [0, e] (e > small enough). By differentiating this identity, we obtain 

—dHixit)) = (graddn (x(t)), x(t)) = (grad„dH{x{t)),x{t)) = 1. 
dt 

Since x = Ph and \Ph \ = I, the eikonal equation implies that ^(^gradjjdn {x{t)),x{t)) = 0, 
or equivalently, that 

gradudn {x{t)) = x(t) 

for every small enough t > 0. This can be written more explicitly as follows: 

Xid„{x{t)) = P, ie Ih = {l,...,h}. (42) 

Let now ct G Iv — {h + 1, ...,?^} and let us differentiate the quantity Xo^dn along the 
normal CC-geodesic x{t) defined by the previous assumptions. We havj^ 

^Xaduix) = {grad(XadH)(x), x) = (gradjj (XadH)(x), x) 
dt 

= {gradu {XadH){x), Ph) = ^ PiXi{XadH){x). 
12 P{0)=Po = P^{0) + P,.{0). 

I'^For sake of simplicity, in these computations we shall drop the dependence on the variable t. 
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Since = + [X„ X„] = X^X, + Y.p^r C^Xp, we get 



X^d^x) = J2PMXdH)ix)+ J2 CtXpd„{x) 
= ^i^,(x„(P,)(x)+ CtXpdu{x)) 



«e/i/ /3e-ford{a) + l 



the first term in the sum is zero because \Ph\ = 1 (unit-speed). Using Definition 11.101 and 
the skew-symmetry of C^, we finally obtain 

j^X^dnix) = - Xpd„{x){CPp„,X^). (43) 

/3G-ford{a) + l 

We summarize the previous discussion in the following: 

Lemma 2.28. Let x : [0, r] — > G (r > 0) he any normal CC-geodesic of unit-speed 
and parameterized by arc-length. Let x{0) = Xq, P{0) = Pq be its initial data and set 
dni^x) = dni^XQ^x) (x G G). Then we have 

(i) gradjidn {x{t)) = Pnit) for every t G [0,r]; 

(ii) igradydnixit)) = -Ef3ei^,,,,,^^,XpdH{x{t))C''P„{t) for every t G [0,r]. 

Proof. The first claim is fj42|) . while the second one is fj43|) . both rewritten using vector 
notation. □ 

We reformulate Lemma 12.281 by using the notation given in Definition 12.51 One has 

^gradydnix) = —C(gradydH(x))PH- 
dt 

At this point Lemma 12.281 can be restated, in geometric terms, as follows: 

Proposition 2.29 (Sub-Riemannian Gauss' Lemma). Let S^{xo,t) the CC-sphere cen- 
tered at xq of radius t G [0,r] and set = T^H\s^(xo,t) a'l^d zu = zu\sn^(^x^^^t) (t G [0,r]). 
Let X : [0,r] — > G (r > 0) be any normal CC-geodesic of unit-speed, parameterized by 
arc-length, with initial data x{0) = Xq, P(0) = Pq, i.e. x(t) = exp s,^{xo, PQ)(t) (t G [0,r]). 
Then, for every t G [0, r] the following O.D.E. 's system holds: 

'^ = -Cu{w)v, (44) 
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Notice that the first equation in (144|) says that each normal CC-geodesic starting 
from xq intersects E>%^{xQ,t) orthogonally (in the horizontal sense), i.e. at the intersection 
point x{t) the velocity vector of the normal CC-geodesic coincides with the horizontal 
unit normal. Furthermore, the second and third equations in (1441) express how change 

= ^H\ssRixo,t) and zu = zJ7|§^(a,o_t) along x{t). 

Proof of Proposition \2.29l We stress that our hypothesis about the absence of abnormal 
minimizers in G implies the smoothness of dn . Since every CC-sphere §^ [xo, t) {t G [0, r]) 
turns out to be defined as E>%^{xo,t) = {x : dn^x) = t}, the Riemannian unit normal 
vector z/ along S%^{xo,t) (at each regular non-characteristic point) may be written just 
by normalizing the following (non unit) normal vector along S^(a;o,t) 

J\f := graddn = {gradudu , gradH^du). 

Clearly, is uniquely determined by A/". The eikonal equation implies that 

— gradHdn, vo = gradydn V t G [0,r]. 

Therefore, by the first equation of system f|T5l) and by (i) of Lemma 12.281 we immediately 
get that x{t) = Puit) = z/^(a;(t)). This identity together with (ii) of Lemma [2.281 implies 
the third equation in (jH]). Finally, the second equation in immediately follows by 
using (i) of 12. 281 together with the third equation of (^^. □ 

Remark 2.30. We stress that Lemma \2. 29\ solves the problem of selecting, for any regular 
point xi belonging to the CC-sphere S^(a;o,'r), the unique normal CC-geodesic having 
velocity vector equals to the horizontal normal direction at that point (i.e. Ph (0) = Uui^xi) ) 
and connecting this point to the center xq of the CC-sphere. Actually, Lemma \2.29[ and 
the uniqueness of solutions of O.D.E. 's imply that the desired curve must be the normal 
CC-geodesic defined by 

x{t) := exp{xi,-Af{xi))(t) t e [0,r] 

where Af = {u^,w). 

Corollary 2.31. Let S = {x & Q : f{x) = 0}, where f is a C'^ function. Assume that 
there exists a CC-ball B{y,r) with center y and radius r, such that B{y,r) C {/(x) < 0}, 
or B{y, r) C {f{x) > 0}, and B{y, r)r\S = {x}, where x ^ S is non- characteristic. Then 
there exists the metric normal 'jj^ to S at x and for every t G [0,r], ■y{t) G 'j^f, where 

7(t) := exps^{y,-^{x)){t) t G [0,r]. 

2.5 2-step case: explicit integration and other features 

In this section we shall explicitly analyze the case of 2-step Carnot groups. Remind that 
in the 2-step setting there exist no abnormal minimizers and our working hypothesis is 
satisfied. In this case, it is well-known that CC-geodesics are smooth; see |34j . 

""^^Remind that, by definition, v^^ := ^| , and m := j^^-^- 
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In order to describe the sub-Riemannian exponential map, we note that the system 
for normal CC-geodesics can explicitly be integrated. Using (fT9|) . i.e. 

X = Ph 
Ph = -Ch{Ph,)Ph 

Ph2 = 0, 

we get that Ph2 G M^^ = i/2 is a constant vector. By standard results about O.D.E.'s, we 
therefore get that 

To obtain the solution in exponential coordinates first note that the equation x = Ph is 
equivalenilll to 

Xi = Pi {i ^ Ih] 

Xa = -\{C'^Xh.,Ph) {aeln^). 

Therefore, setting XH{t) := x/i(t)} G to denote the projection of the 

solution x{t) of (|T9l) onto the firsts h variableqlj, one gets 



Xh it) = Xh (0) + [ e-^" (^^2 (0) ds 
Jo 



(45) 



and 



1 /■* 

Xa{t)=x^{0)-- {C'^x„,XH)ds. (46) 
^ Jo 

These equations describe the sub-Riemannian exponential map in the 2-step case. 
More precisely, fixing a base point Xq, we have that 

exp^(xo, ■)(■) : UH X H2XR — >G 

is given by 

exp^(a;o,Po)(t) :=xo + ^ e-^«(^«2)^PH(0)ds-^ | ^ (C^Xh, Xh) rfsje,. (47) 



aelH2 

Here above UH denotes the bundle of all unit horizontal vectors, i.e. UH (Z H = E>' 



h-l 



-"^^Explicitly, we have 



'^^Using exponential coordinates for G, every point a; G G is n-tuple x = (xi, x^, Xh+i, Xn). So it 
seems natural to "divide" the variables as follows: 

Xh ■■= (xi, ...,Xh) e M'\ XH2 := {xh+1, ■■■,Xn) G M^^ (n = /l + /12) 
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Remark 2.32. Setting x(t) := exp gj^{xo, Po)(t) and 

£{t) := f e-^^^^'^'^^'ds e Mh, 
Jo 

we get that 

x{t)=xo+l^{S{t)VH{0),ei),...,{S{t)VH{0), en), (C^i^, x„) rfs) , . . . ) , 



a— th place 

Notation 2.33. In the sequel, we shall denote by dxn the following vector valued 1-form: 

dxH ■■= {dxi, dxh)^' G H* x . . . x H*^ . 

h — times 

Notation 2.34. Let x : [0, T] — > G be a CC-normal geodesic such that x{0) = xo and 
x{T) = xi. Later on we shall set [xo,Xi] := s([0,T]) = {x{t) G G : t G [0,T]}. 

Remark 2.35. For any a G 1^2 ; the integral 



{C ^ X H •, X ds 

'[xh{0),xh (*)] 

can explicitly be evaluated, by means of standard linear algebra arguments, by noting that, 
as every skew- symmetric linear operator, can be written, after an orthogonal change 
of basis, in a "canonical" form. More precisely, there exists G Oft(]R)0 such that 

/ A° ... \ 

-A? 











-A^ 



A^ 



V 



where ±iA" (j = are the eigenvalues - purely imaginary - of computed with 

their multiplicity, and On denotes the zero N x N -matrix where N is the nullity of ■ 
Here R = |rankC" and h = 2R + N. Setting 

ynit) = (0")-ix4t), yniO) = (O")-1xh(0), 



we get that 
it) 



R 



{C^Xh , dXn ) 



XH(0),XH{t)] 



[VH {0),yH (t)] 



yj+idyj - yjdyj+i. 



We stress that similar integrals was previously introduced and studied by Pansu in his deep 
study of differentiability in CC-spaces; see f5^ . Definition 4-6, p. 15. 



^^i.e. the orthogonal group on M''(= H). 
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Using a standard result about periodic solutions of linear O.D.E.'s0, we may study 
the existence of T-periodic solutions, T > 0, of the h-th order linear O.D.E.'s system 



Ph = -Ch{Ph,)Ph, 



(49) 



where Ph2 G H2 is constant. More precisely, we get that there exists a /c-dimensional 
sub-space Sk of the solution space Sh of the equation (l49l) such that each member of Sk 
is T-periodic, and no member of Sh \ Sk is T-periodic, if and only if the following holds: 

rank(e~^^ ^^"^ - Id/,) = h - k. 

Hereafter, we shall analyze this condition. First, note that Ph2 = implies that Ph 
is constant and so we may assume Ph2 7^ 0. As in Remark 12.351 we make use of a 
standard Linear Algebra argument. More precisely, as every skew-symmetric linear oper- 
ator, Ch {Ph2 ) can be written, after an orthogonal change of basis, in canonical form. In 
particular, there exists O G 0/i(]R) such that 

\ 



Ch{Ph2) ■.= 0-^Ch{Ph2)0 






Ai 








-Ai 




















A2 








-A2 














V 0^ / 

where ± i Aj (j = 1, i?) are the - purely imaginary - eigenvalues of Ch (-Phj ), computed 
with their multiplicity, and On denotes the zero x A^-matrix, where A^ is the nullity 
of Ch{Ph2). Here R = |rankC// (Pi/2 ) and h = 2R + N. Notice that the eigenvalues are 
functions of Ph2, i.e. Xj = Xj{PH2) (j = 1,...,R). Setting 

Pn{t) = 0-'PH{t), PniO) = 0-'PHiO), 
we obtain the following equivalent equation: 

Pn = -Ch{Ph2)Ph. 



(50) 



^^The following result holds true (see [5]): 
Theorem 2.36. Let 



ij := A{t)y 



(48) 



be a n-th order system of linear homogeneous differential equations, where A{t) = [aij(t)]ij=i,...,Ti is a 
n X n matrix and each element aijit) is a real-valued function continuous on the real line M. Let W 
denote a fundamental matrix of (I48|) . Let k he a non-negative integer, < k < n. Then there exists 
a k-dimensional sub-space Sk of the solution space Sn of (|48|) such that each member of Sk is periodic 
of period T and no member of Sn \ Sk is periodic of period T if and only if the rank of the matrix 
W{T) - W(0) isn-k 
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By applying the previous argument, we infer that there exists a /c-dimensional sub-space 
Sk of the solution space Sh of the equation (!50!) such that each member of Sk is T-periodic, 
and no member of Sh \ Sk is T-periodic, if and only if: 

rank(e-^ ^^"^ - Id/,) = h - k. (51) 

Therefore, it remains us to analyze the matrix e~^" ^^"^ ^'^ given by 



/ cos(Ai(Ph2)T) sin(Ai(PH2)r) 0... 
- sin(Ai (Ph2)T) cosI\iIph^)T) 0... 



: ; cosi\RiPH,)T) smi\RiPH,)T) 

-sm{XR{PH,)T) cos(Afi(PH,)T) 
V Id^ / 

We claim that for every given Phj 7^ 0, there exist T-periodic solutions of (H9l) . for 
some positive T = T(Ph2). Indeed, if 7^ 0, i.e. the nullity of Ch{Ph2) is non-zero, 
the claim follows by the above discussion since rank(e~'"^ ^^"2 — id^) > h — N and 
k > N > 0. Furthermore, if = 0, the claim follows by choosing T = 2tt / Xj{PH2) 
for some j = 1,...,R. Indeed in such a case one sees that k > 2. With an analogous 
argument, one can easily show that the dimension of the space Sk of T-periodic solutions 
of (Sni) satisfies: 

N + 2<k<h. 

Note that these solutions can also be constant functions. In similar way, we see that for 
every Phj 7^ there exists T = T(Ph2) > and there exist -at least- two non-constant, 
linearly independent, T-periodic solutions of ( l49i) . 
We summarize the above discussion in the next: 

Proposition 2.37. For every Ph2 G H2, Phj 7^ 0, there exists T = T(Ph2) > and 
there exists a k-dimensional sub-space Sk of the solution space Sh of (1501) such that each 
member of Sk is T-periodic and no member ofSh\Sk is T-periodic. Furthermore it turns 
out that N -\-2 < k < h, where N is the nullity of Ch{Ph2). Finally, there exists a positive 
integer r > 2 and there exists a r -dimensional sub-space Sr of Sk such that each member 
of Sr is T-periodic and non-constant. 



Remark 2.38. In the previous Proposition\2. 37\ the numbers k and r can be characterized 



in a slightly different way. Indeed it turns out that k is the multiplicity of the eigenvalue 
X = 1 of the scalar matrix e~'"'^*-^^2)'^. Moreover r = k — N . To prove this claim one 
can use another characterization of T-periodic solutions of linear homogeneous systems 
of O.D.E. 's which can be found in 19]. 

We may use the above Proposition 12.371 to study the T-periodicity of Xh] see (H5ll . 
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Remark 2.39. If Xh is T-periodic then PH{s)ds = 0. This follows by hypothesis, 
using the first equation of (fT9|) . Note also that if Ph is T-periodic, then Xh is T-periodic 
if and only if Ph (s) ds = C0. 

Example 2.40. Let Ph2 ^ and let us assume that Ch{Ph2) is invertible. For instance, 
this is the case of the Heisenberg group H" . In such case, by a direct computation based 
on the very definition o/e"'"^*-^"^)'^^ a turns out that 

Xh (t) = Xh (0) + {Ch {Ph, (id, - e"^" ^^^^ )*)p^ (Q). 

Remark 2.41. Let us consider the system of normal CC-geodesics (|T9l) . We would like 
to remark that the T -periodicity of Xh is related with some other things. To this aim, let 
us assume that Xh be T-periodic. More precisely, we assume that there exists a minimal 
T = T{Ph2) > such that Xnit) is T-periodic, for any given Ph, ^ 0. Furthermore, let 
xi := exp sn.{xQ, Pq){T) , where Pq = (Ph (0), Pf/j ) G UH x H2. Then it can be shown that: 

(i) The point xi is conjugate to Xq along the normal CC-geodesicx{t) = exp g,^{xo, Po){t), 
te[0,T]; 

(a) The "segment" [xo,Xi] := a;([0,T]) = {x(t) G G : t G [0,T]} is a minimizing normal 
CC-geodesic and Xi is the cut-point of Xq along the normal CC-geodesic x{t). 

In this case, we will say that [xo,Xi] is a minimizing CC-geodesic segment. Note 
that the CC-length of segment [xq, Xi] is simply T . The proofs of these claims can be done 
by following a classical pattern, for which we refer the reader to JW^ . However, this is 
beyond the scope of this paper. 

Remark 2.42. For 2-step Carnot groups, similar arguments can be used to show that the 
CC-sphere §^(a:;,r) = {y E G : dH{x,y) = r} is -smooth out of the set 

H2{x)nSl,{x,r), 

where 

H2{x) := {z = exp^{zH , zh,) e G : Zh = Xh}- 

More precisely, each point y G (x, r) \ H2{x) can be joined to the center x of SJ5^ (x, r) 
by a unique minimizing normal CC-geodesic. In particular, one can show that for every 
y G S^(x, r) \ H2{x) there exists a unique Pq G UH x H2 such that y = exp g,^{x, Po){r) 
and dexp sTi{x, Pq) has maximal rank. On the other hand, if y & §5i^(x, r) fl H2{x), then 
it turns out that y G Conj(x). 

2.6 Appendix: iterative integration for normal CC-geodesics 

In this section we will how how, at least in principle, the system of normal CC-geodesics 
canbe integrated, step by step. In what follows, we shall use the notation Ph^ = Vh^ (P). 
Remind that (JTSll is given by 

j X = Ph ( <^==^ Xh = Ph, Xv = O) 

\ P = _C{Pv)Ph] 

^^Remind that if / : K — > R is a continuous T-periodic function, then /(s) ds is T-periodic if and 
only if f^ f{s) ds = 0. 
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By orthogonal projection onto Hk, we get that Ph,^ = 0. Hence Ph^. is constant. Then the 
k — 1-th vector equation becomes Ph^^^ = —C{Ph,.)x and so 

Ph,_, (t) = Ph,_, (0) - C{PH,){x{t) - xo). 

By iterating the same procedure one gets Ph^_2 = —C{P h,,^^)xh ■ Hence 



P».-2 {t) 
Ph2 {t) 



Ph,.2{0)- [ C{PH,_,{s))Xn ds 

Jo 

Ph,{0) - [ C{Ph,_,{s))xh ds 
Jo 



Ph2(0)- / C{Ph,{s))xh ds 
Jo 

Ph(0)- / C„{PH,{s))xHds. 



Finally we get that 

XH{t) 



Xh(0)+ f PH{s)ds 
Jo 

Xh (0) + ^ (Ph (0) - _^ ' C'h (Ph, {s2))xh ds2^ ds, 



(0) + Ph (0)t 



t fSl 



^ JO 



Ch (Pff2 {s2))xh ds2 ) ds 



or equivalently, that 



Xh {t)=CH{PH,{t))XH. 



In the sequel we shall apply this procedure to the case of general 3-step Carnot groups. 
We stress that the first vector equation in (I2T!) . i.e. a: = Ph, in exponential coordinates, 
is equivalent!^ to the following system: 



X. 



_ P. 



[i G Ih) 

< Xa = -\{C'^Xh.,Ph) (aG/Ha) 

^ ip = -{\{C^x,Pu)-^^T.o.^j^ {C^x,e^){C^Xn,PH)) iPeln,] 



^"By the formula for left-invariant vector field stated in Example 1 1.24i we get that 



ieln 



ieln 



Xi3 



{a G Ih2 ) 
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Therefore, for any a G /hj , we have 



Xa{t) = Xa{0) - - 



{C^Xh , dXn ) 



(52) 



XH {0),XH (t)] 



and, for any P E Ih^, we have 

Xfj{t) = Xf3{0)- (^^{C^x,Xh) - ^ ^ {CHXH,XH){C'^x,ea)^ds 

:,dxH)-^ {C^x,ea){CHXH ,dxH)^ (53) 



Xf^iO) 



xh{0),xh (t)] 



Now we proceed with the integration of (pT!) as we have explained in general, at the 
beginning of this section. As before, we have that Ph^ = and so Ph.j is constant. Since 

= -C{Ph.^)xh we get 

PH,{t) = Ph,{0) - C{PH,){xHit) - Xh{0)). 

Moreover Ph = —C{Ph2)xh, and so 

Pnit) = Ph{0)- [ CH{PH,{s))xH{s)ds 
Jo 

= Ph(0)-/ CH{PH,iO)-C{PH,)[xHis)-XHiO)])xHis)ds 



Finally, we have 



Xh (t) = Xh(0)+ / PHis)ds 



XHiO)+ I (P«(0) 



= Xh{0)- 

= Xh(0) +PH(0)t- 

or, equivalently. 



Ch {Ph2 {s2))xh (S2) ds2j dsi 

Ch {Ph,{0) - C{Ph,)[Xh{s2) - Xh{0)])xh{s2) dS2) dsi 



Ch ( Ph2 (0) — C(Ph3 ) [xh (S2) — Xh iO)])xH ds2 



dsi 



X 



^ = -Ch{Ph2{0) - C{PH,)[xHit)-XHiO)])xH. 



(54) 



=Ph2 it) 
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3 CC-distance from hyper surfaces in 2-step Carnot groups 



Throughout this section we shall assume that G be a 2-step Carnot group. Let S" C G be 
a smooth hypersurface (i.e. closed (n — l)-dimensional submanifold of G) and let 



where is the Riemannian unit normal along S. We stress that the normal (non-unit) 
vector field A/" is defined at each non-characteristic point x & S \ Cs- 

Remark 3.1. We would like to remind some classical results about the regularity of the 
distance function to smooth hypersurfaces (or, more generally, submanifolds) which hold 
in the Euclidean -and Riemannian- setting. If N G is a C^-smooth manifold with 
k > 2, then one easily sees that, near N, the distance function 6^ is C''^^ -smooth. A 
first, somewhat surprising result, was proved by Gilbarg Trudinger in the Appendix of their 
celebrated book I3^ . Indeed they prove the existence of a neighborhood U of N such that 
the distance function Sn is of class C*^ on U \ N. A similar problem was also considered 
by Federer in its theory of Sets of Positive Reach; see f2^ . The complete solution of this 
problem, in the case, can be found in a paper by Krantz and Parks; see l4^ . We stress 
that, in such a case, a further hypothesis is needed. More precisely, if we assume that N is 
just -smooth and that there exists a a neighborhood U of N having the so-called "unique 
nearest point property" , then 5^ is of class on U \ N. Later on, a simplified proof of 
this result was proved by Foote; see 124]. We would like to remind that, some previous 
results of this type, in a sub- Riemannian setting, was proved by Arcozzi and Ferrari for 
the first Heisenberg group M.^ ; seef^. 

So let us begin by assuming that S is C'^-smooth with k > 2. As a first thing, let us 
define the mapping $ : Sx] — e, e[ — > G (e > 0) 



By construction, it turns out that $ G C''^^{S \ — ^,^[) for any - small enough - 

e > 0. The "size" of e depends on the local geometry of S. 

Next we shall apply the Inverse Mapping Theorem to this map. For this reason, we 
have to compute the differential of $ at t = 0. To this aim, for any fixed point yo G S, 
let us introduce a system of Riemannian normal coordinates^ in a neighborhood of yo; 

^^More precisely, e represents the CC-length of a minimizing CC-geodesic segment having "lagrangian 
multipher" Ph2 = '^{y)\ see Remark 12.411 

^^Let (ri, ...,T„_i) be an orthonormal basis of Ty^S. For 1 < j < n — 1 define a real-valued function 
Uj on a neighborhood of yo by 



where exp^ denotes the Riemannian exponential map. Then, by definition, (ui, ...,u„_i) is a system of 
normal coordinates corresponding to the orthonormal basis (ti, ...,r„_i). 



5h -.G — ^ U {0} 

denote the CC-distance function for S, i.e. 6h (x) = infyg^ dn {x, y). 
In the sequel we shall use the notation 




$(l/,t) :=exp^(y, Afiy)m. 



(55) 
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see, for instance, [TT] or [39]. Hence y = y{ui, and 

d d 

dui ' ' dUn-l 

are the coordinate vector fields associated with the normal coordinate system (ui, Un-i)- 
Using these coordinates, it turns out that V{y) := ^ A ... A ^ is a normal (non-unit) 
vector along S*, in a neighborhood of j/o G S. By definition, the Riemannian unit normal 
v along 5* (in this neighborhood of yo) is given, up the the sign, hy v = -j^, while the 

horizontal unit normal z/^ can be written out as z/^ = ; see Section 11.21 

Lemma 3.2. Let S C G be C^-smooth with k > 2. Then 

\det[j(^y,o)'^]\ = \VHV\, (56) 

where J{yfi)^ denotes the Jacohian matrix operator at {y,0) G Sx] — e, e[. Therefore 
I det [j7(j/,o)'^'] I 7^ at each non- characteristic point of y G S \Cs- 

Proof. Throughout this proof we will need some results of Section 12.51 By ( HTIl we get 



^y,t)=y+ I e-^-(-(^))V,(y)rf.-l (/ 



{C^Xh (s), Xh (s)) ) e^. 



where 



Jo 

+ / e-^-(-(^))V,(i/)rf.),e.), 



To sake of simplicity, we also set 

^y,t) ■.= y + Aiy,t). 

By using the explicit expression of $ and the Fundamental Theorem of Calculus we ge^^'^' 

dt 



dt 



- {C^ X H , X h) Qa- 
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Since 2:(i) := ^{y,t) G G, in exponential coordinates we have $ = expg(($/f , $/f2 )) = (<i>if , $if2 ), 

if2 = $i/2 



where = G M'' and xh2 = <f>i/2 G M^'; see Notation [Tl 
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In the last formula we have also used the skew-symmetry of the structure constants 
operators (a G /hj ) . Note that from the previous computations we obtainQ 



t=0 



Tr 



= ^niy)- 

Furthermore, using the system of normal coordinates u = {ui, ...,m„_i), we get that 



a$ _ dy OA 

dui dui dui 



As usual we will set A = exp MAh , AH2)) = {Ah,Ah2). We compute 



dAn 

dui 

OAh^ 
dui 



OUi 



grady{{CHXH , Xh)) , ds 



{i = 1, - 1). 



Therefore, choosing t = 0, one gets 

J{yfl)^ = col 

= col 



dy 9$ 



dy 

dui ' ' dun-i ' dt 
dy dy 



dur 



t=OJ 



VdUi uun-l 

Finally, we may compute the Jacobian determinant of Ji^y^-^^. By using standard Linear 
Algebra arguments, one gets 

dy dy 



det [J{yfi)^] 



det l^col 

dy 



dui ' ' dUn-l ' " 

dy 



dui 



A ... A 



dy 



dui 
dy 



A ... A 



A ... A 



dUn-l 

dy 



dUn-l 

dy 



dui dUn-l 

VnViy)], 



My),'^Ay))\ 

\Vnu{y)\ 



which achieves the proof. 



□ 
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Remind that we are working in exponential coordinates. So we have 

is/if iS/ff a£lHo 
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Remark 3.3 (Invertibility at the non-characteristic set). Let us set So '■= S \ Cs- It 

turns out that Sq is an open subset of S , in the relative topology. Moreover, since we are 
assuming that S is -smooth with k > 2, one gets that dimC5 < (n — 2); see \4^ - 
Now let Uq (s 5*0 be an open set compactly contained in Sq. By Lemma \3T^ we know that 
the Jacobian of the mapping $ : ?7ox] — e, e[ — > G is non-zero along Uq x {0}. The Inverse 
Mapping Theorem implies that there exists g]0, e] such that 

$ : f/ox] -eo,eo[ — ^ <l>(f/ox] - eo, eo[) 

is a C'^^^ -diffeomorphism. 

Notation 3.4 (Projection mapping). The previous Remark \ 3. 3\ enables us to define the 
following mapping: 

^ := $-1 : ^{Uox] - eo, eo[) — > Uox] - eo, eo[. 

By construction, is Q^^^-smooth. In the sequel, we shall denote by ^ s the projection 
of the map ^ onto its 1st factor, i.e. \E'(a;) = (\[^5(x), t(x)). 

Let us set U := ^{Uqx] — eo,eo[) C G and let x E U. The previous discussion can 
be summarized by saying that every open set Uq, which is compactly contained in 5*0, 
has a neighborhood U C G satisfying the unique nearest point property^^ i.e. for every 
X eU there exists a unique point y E Uq d Sq such that 5h{x) = dn^x^y). By using the 
previous notation, one has \l/(x) = (y,t), where y = \E'5(x) and t{x) = dn^x^y) = 6h{x). 

Theorem 3.5. Let G be a 2-step Carnot group. Let S G G be a C^-smooth hypersurface 
with k > 2 and let 6h denote the CC-distance function for S , i.e. 6h (x) = infyg5 dn {x, y). 
Set Sq := S\Cs, where Cs denote the characteristic set of S. Then, for every open set Uq 
compactly contained in Sq, there exists a neighborhood U <Z G of Uq having the unique 
nearest point property with respect to the CC-distance dn ■ Finally, the CC-distance 
function from UqH S is 6h \ u\Uo C'' -smooth function. 

Proof. We just have to prove the last claim. To this aim, let X G X(G) and set y := \I's'(x), 
where ^^(x) denotes the projection along Uq C 5*0 of the point x eU := ^{Uqx] — eo, eo[). 
Moreover set t := dni^x^y). We have 

{graddn (x), X) 
= {graddnix,'^ s{x))-,X) 

= {{9radJ„{x,y))\^^^^^^yX) + (^[j.^'^s{x)]X, {gradyd„ {x,y))\y^^^^^^y (57) 

Now let us introduce the following further notation. 

Remark 3.6. Let x E G, y E G \ H2{x) = {z = exp^^Zn , ZH2) E G : Zh = Xh} and set 
t := dH{x,y). Moreover, let 

^^With respect to the CC-distance dn . 
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denote the Riemannian unit normal along the CC-sphere Ssiz {x, t) at any regular point 
y G Ssiz{x,t). We stress that, under our hypotheses, the point y turns out to be a regular 
point of Eisjz{x,t) ; see, for more details, Remark \2.4^ Below we shall set 



According to the results of Section [2^ we immediately get that: 

grad^daix.y) = A/'y,t(x), 
gradydH{x,y) = N'x,t{y)- 

By applying the results of Section 12.41 together with the explicit form of CC-geodesic for 
2-step Carnot groups (see Section [275]) . we get that 

Ar,,(x) = (e-^-(-(^))V,(i/),ti7(y)). 

Notice that My^t is C^~"^-smooth as well as A/" = (z/^,ti7). Furthermore, it is easy to see 
that Mx,t{y) = =tA/"(?/), where the sign only depends on the given orientation of S. 
By the previous discussion and ( 1571) we get that 



{gradSH{x),X) = {Af^^(^,^^t{x), X) + {[j,^ s{x)]X,Af,^t{^ s{x))). 

Since [j72;\E's(x)]X G T^g(^x)S, by using the fact that M^^ti"^ si^)) is normal to S at \E'5(x) 
one gets 

([Xvl/^(x)]X,A4,,(vI/s(a;))) =0. 

Therefore {gradSn {x),X) = (A/*<;(x'),t(x), X). By the arbitrariness of X G X(G), it follows 
that gradSn is of class C''^^ on U \ Uq. Hence^H is of class C'^ on f/ \ Uq. This achieves 
the proof. □ 
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